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ABSTRACT 


Because  a turbulent  boundary  layer  in  a nozzle  undergoes  strong 
acceleration,  a 1 ami nari scent  boundary  layer  with  the  benefit  of  a 
lower  surface  heat  transfer  rate  can  result  for  some  short  distance. 
Here  and  in  a previous  report  (Simpson  arid  Wallace,  1975)  several 
phenomena  which  accompany  laminariscence  produced  by  strong  acceler- 
ation are  examined  for  nozzle- flow  and  sink-flow  accelerational  dis- 
tributions, respectively.  Several  uifferent  type  measurements  of 
the  structure  of  two  nozzle-type  flows  are  reported  to  determine  how 
an  initially  normal  turbulent  boundary  layer  approaches  the  laminar- 
like  state,  including  mean  velocity  and  Reynolds  stresses  profiles, 
spectra,  turbulent/non-turbulent  interfacial  structure,  and  wall 
bursting  and  sublayer  spanwise  spatial  structure. 

As  a result  of  these  experiments,  it  appears  that  the  surface 
skin-friction  is  not  reduced  to  laminar  values  in  sink  flows  unless 
K(  =vUoo'cl l^/dx)  is  greater  than  about  3.6  x 10'6.  In  nozzle-type 
flows?  K°°must  also  be  greater  than  this  value  over  a short  distance 
in  order  to  produce  a short  rel ami nari zed  region  downstream.  The 
large-eddy  structure  of  the  outer  region  governs  the  bursting  fre- 
quency, the  intermittent  bulge  passage  frequency,  and  influences 
the  wall  flow  behavior  downstream.  After  the  cessation  of  entrain- 
ment of  free-stream  fluid,  these  frequencies  approach  constant  values. 
The  wall  spanwise  structure  appears  to  lag  behind  local  conditions 
and  to  reflect  upstream  flow  behavior.  After  retransition  to  a low 
acceleration  turbulent  boundary  layer  downstream,  much  larger  spanwise 
scale  structures  are  observed. 

The  entrainment  rate  of  non-turbulent  fluid  decreases  to  zero 
at  about  the  streamwise  location  at  which  the  shape  factor  reaches 
a minimum  value.  The  cessation  of  entrainment  by  the  eruption  and 
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and  engulfment  action  of  the  large  eddies  can  be  traced  to  the  reduc- 
tion of  available  turbulence  energy  for  diffusion  to  free-stream 
fluid.  This  reduction  of  available  energy  is  due  to  the  negative 
normal  stresses  turbulence  energy  production  term.  Spectral  distri- 
butions of  the  streamwise  fluctuation  F(n)  possess  a frequency  region 
where  nF(n)  is  constant  for  laminar-like  boundary  layers  at  large 
K values. 
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NOMENCLATURE 


A,  B,  C,  D 
Cf/2 
E,  e 
F(n) 

F’  F6 
f,  f' 

H 

K 


constants. 

2 

= T^pU^,  friction  factor, 
mean  and  fluctuation  voltages, 
spectral  density  defined  by  equation  (9). 
defined  by  equations  (18)  and  (20),  respectively. 
U/U^  and  9f/9n»  respectively. 

= 6*/0,  shape  factor. 

dUoo 

= ; also  a constant  in  equation  (12). 

^oo 


M 

m 

n 


ni 

dP/dx 

«i 

7 

sR  (z,  T ) 

XT ' S' 


sR(z) 

Ree 

Re. 


integral  length  scale  defined  by  equation  (14). 
constant. 

exponent  in  split-film  calibration  curve, 
frequency,  Hz. 

frequency  of  intermittent  bulge  passage  at  Y. 
pressure  gradient. 

power  dissipation  in  split-film  sensor  i. 

~2  ~2  ~Z 

= u + v + w . 

normalized  spatial  cross-correlation  of  two  x surface 
fluctuation  signals  during  sample  time  T . 

A 

sR(z,  Ts)  as  Ts  becomes  very  large. 

= U^e/v,  momentum  thickness  Reynolds  number. 


, microscale  Reynolds  number. 


length  of  sample  time, 
time  delay,  time. 


' 


i 


1 


U,  V,  w, 

U,  V,  W 

-u7 

VE 

x,  y,  z 


Y 

SUBSCRIPTS 

b 

eff 

i 

a 

w 


■ (tw/p)^2,  shear  velocity. 

velocity  fluctuations  in  the  streamwise,  normal,  and 
spanwise  directions. 

mean  velocities  in  the  streamwise,  normal,  and  span- 
wise  directions. 

kinematic  Reynolds  shearing  stress, 
entrainment  velocity. 

cartesian  co-ordinates  in  the  streamwise,  normal,  and 
spanwise  directions,  respectively. 

distance  from  the  wall  to  where  y = 0.5. 

denotes  "bursting"  value, 
denotes  effective  cooling  velocity, 
split-film  sensor  index, 
denotes  linearized  signal, 
denotes  wall  value, 
free-stream  condition. 


y 

a 


GREEK 

a 


Y 

6 


constant  in  equation  (10);  split-film  probe  yaw 
parameter. 


intermittency,  long-time  averaged  fraction  of  time 
that  the  flow  is  turbulent. 


60  99  = y where  f = 0.99;  60;9g5  = y where  f = 0.995. 
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dy,  displacement  thickness. 


dy,  momentum  thickness. 
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n 


— 5"  1 /2 

= (15vu  /e)  ' , a turbulence  microscale. 

spanwise  wavelength  in  wall  sensors  cross-correlations. 

defined  by  equations  (15). 

kinematic  viscosity. 

density. 

standard  deviation  for  intermittency  distribution, 
shearing  stress. 

angle  of  flow  incidence  to  split-film  probe, 
logarithmic  equation  (5)  constant. 
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The  strong  acceleration  of  turbulent  boundary  layers  in  a nozzle 
can  produce  the  benefit  of  a substantially  reduced  convective  heat 
transfer  rate  when  relaminarization  occurs  (Moretti  and  Kays,  1965). 

The  accelerational  parameter  K,  defined  as  vU^  dU^/dx,  must  increase 
to  a value  in  excess  of  about  3.6  x 10~^  (Kline  et  al . , 1967)  for 
relaminarization  to  occur.  The  practical  consideration  that  flow 
acceleration  extends  over  a finite  length  then  requires  that  K decrease 
to  zero  downstream.  Narasimha  and  Sreenivasan  (1973)  reexamined  the 
behavior  of  this  type  of  flow  from  the  time-averaged  results  of  many 
previous  experiments.  However,  there  is  little  information  on  the 
flow  structure  that  would  explain  the  mechanisms  for  relaminarization. 
Exceptions  are  the  space-time  correlations  of  Blackwelder  and  Kovasznay 
(1972)  for  the  large-scaled  motion  and  the  wall  region  structure 
reported  by  Schraub  and  Kline  (Kline,  et  al . , 1967). 

The  work  described  here  is  the  second  part  of  a program  to  pro- 
vide experimental  flow  structure  information  for  strongly-accelerated 
turbulent  boundary  layers  under  different  K distributions.  In  the 
first  part  (Simpson  and  Wallace,  1975)  results  were  obtained  for  two 
asymptotic  sink-flow  type  turbulent  boundary  layers  with  constant  K 
distributions  (flow  A,  K = 2.17  x 10-6;  flow  B,  K = 3.19  x 10"6).  Here 
results  are  presented  for  a flow  in  which  K increases  linearly  with  x 
(flow  C)  and  a flow  in  which  K increases  linearly  with  x to  a maximum 
value,  stays  constant  at  this  level  for  a short  distance,  and  then 
decreases  linearly  with  x (flow  D).  Flow  C models  the  rapid  increase 
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in  K encountered  in  nozzles  while  flow  0 models  both  the  increase  and 
decrease  in  K.  Since  the  acceleration  of  both  flows  took  place  in  the 
same  length  of  wind  tunnel  test  section,  dK/dx  for  the  increasing  K 
portion  of  flow  D was  more  than  double  dK/dx  for  flow  C. 

The  type  of  measurements  reported  here  for  flows  C and  D are  simi- 
lar to  those  obtained  for  flows  A and  B.  In  addition  to  mean  velocity 
and  Reynolds'  stresses  measurements,  spectra  of  the  streamwise  fluc- 
tuations were  obtained.  This  was  partly  motivated  by  the  presence  of 
a flat  nF(n)  spectrum  for  flow  B and  the  question  as  to  whether  this 
behavior  accompanies  true  relaminarization.  Other  measurements  are 
the  turbulent/non- turbulent  interfacial  intermi ttency  and  frequency 
of  passage  of  intermittent  turbulent  bulges,  the  wall  "bursting" 
frequencies  or  rate  of  passage  of  eddies  over  the  wall,  and  the  span- 
wise  spatial  structure  at  the  wall.  This  latter  group  of  measurements 
provides  further  insight  as  to  the  roles  that  the  bursting  behavior, 
spanwise  structure,  and  intermittency  play  in  relaminarization.  Since 
these  measurements  have  been  obtained  for  four  different  K distribu- 
tion flows  on  the  same  apparatus,  the  effect  of  the  K distribution  on 
the  developing  flow  structure  can  be  determined. 

2.  EXPERIMENTAL  APPARATUS 

In  general,  all  of  the  apparatus  and  instrumentation  described 
by  Simpson  and  Wallace  (1975)  for  flows  A and  B were  used  for  flows  C 
and  D.  The  SMU  wind  tunnel  with  a sixteen  feet  long,  three  feet  wide, 
test  section  was  used  to  produce  the  desired  boundary  layer  on  the 
flat  bottom  wall  by  adjusting  the  plexiglas  top  wall.  In  the  current 


experiments  the  freestream  flow  upstream  of  the  acceleration  was  uni- 
form within  0.066%  in  the  spanwise  direction  and  within  1%  in  the 
vertical  direction,  with  a streamwise  turbulence  intensity  of  0.6%  at 
9.1  fps.  At  the  tunnel  exit  the  free-stream  streamwise  turbulence 
intensity  dropped  to  below  0.55%  and  0.5%  for  flows  C and  D,  respec- 
tively. Of  this  intensity,  about  0.5%  was  due  to  unsteadiness  at  about 
20  Hz.  Figure  1 is  a side  view  schematic  of  the  test  section  with 
the  upper  wall  locations  for  flows  A,  C,  and  D.  In  all  four  flows 
the  flat  upper  wall  was  15  inches  above  the  test  wall  at  the  entrance 
and  19.5  inches  above  it  at  the  96  inches  loca  ion.  For  flow  C the 
parabolic-shaped  upper  wall  sectior  began  at  9b  inches  and  was  1.72 
inches  above  the  test  wall  at  the  exit  at  1 90. 8 inches.  For  flow  D 
the  antisymmetric  ogee-shaped  upper  wall  began  at  96  inches  and  was 
1.84  inches  above  the  test  wall  at  189.7  inches.  All  experimental 
data  were  obtained  with  the  temperature  maintained  constant  at  77 
±1/2°F  and  a constant  stagnation  pressure  at  the  exit. 

The  1/4  inches  blunt  leading  edge  trip  on  the  test  wall  and  the 
boundary  layer  smoke  injection  arrangement  upstream  of  the  trip  were 
the  same  as  used  for  flows  A and  B.  The  same  smoke  generation  system 
with  a mean  particle  size  of  about  1 micron  was  used  when  smoke  was 
used  to  mark  turbulent  fluid  upstream  of  relaminarization.  For  these 
intermi ttency  measurements  the  optics  and  traversing  equipment  of  the 
SMU  laser  anemometer  were  used  as  for  flows  A and  B. 

Standard  Thermo-Systems,  Inc.,  model  1050  constant  temperature 
anemometers,  model  1055  linearizers,  model  1057  signal  conditioners, 
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and  model  1015C  correlator  were  used.  A TSI  model  1274-10  boundary 
layer  hot-film  probe  was  used  in  measuring  the  mean  and  streamwise 
fluctuation  velocities  in  flow  C while  a model  1218-T1.5  boundary 
layer  hot-wire  probe  was  used  for  flow  D.  The  sensing  element  for 
the  hot-film  is  a 0.001  inches  diameter  platinum  coated  quartz  rod 
with  a sensing  length  of  0.04  inches.  The  0.00015  inches  diameter 
platinum-plated  tungsten  wire  had  a 0.05  inches  sensing  length.  Based 
on  the  Col  1 i s and  Williams  (1959)  equation  for  the  anemometer  bridge 
output  calibration,  the  linearizers  were  adjusted.  The  linearizer 
output  was  directly  checked  with  the  known  calibrator  velocity  for 
accuracy  and  linearity,  the  maximum  tolerable  deviation  being  less 
than  about  1/2%.  Consequently,  the  uncertainty  on  velocity  measure- 
ments is  about  ±1/2%. 

The  split-film  probe  (TSI  model  1287)  was  used  to  determine  U, 

T ~Z  

V,  u , v , and  -uv  for  flow  C.  This  relatively  new  hot-film  probe 
was  selected  because  of  its  relatively  small  size  since  the  boundary 
layers  under  study  were  thin.  The  split-film  sensor  is  a modifica- 
tion of  the  basic  platinum  coated  cylindrical  film  sensor.  Two  elec- 
trically independent  films  each  cover  one-half  of  the  circumference 
of  a 0.006  inches  diameter  quartz  rod.  Each  film  is  operated  by  a 
separate  constant  temperature  circuit.  The  non-uniform  heat  transfer 
distribution  around  a constant  temperature  cylinder  is  utilized  to 
measure  the  fluctuating  components  of  the  instantaneous  velocity  vec- 
tor. To  avoid  thermal  variations  in  the  substrate  the  sensors  must 
be  held  at  closely  identical  temperatures.  The  resulting  output 
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voltages  are  used  in  the  same  manner  as  those  from  an  x-wire  probe  to 
determine  mean  and  fluctuation  quantities. 

Using  a right-handed  coordinate  system,  the  plane  of  the  two 
splits  that  separate  the  two  platinum  films  from  one  another  are  in 
the  x-z  plane,  being  nominally  parallel  to  the  test  wall  in  these 
experiments.  Following  Spencer  and  Jones  (1971),  the  power  dissipated 
in  each  film  can  be  related  to  the  velocity  by 

Q,  - (A  * Bu”ff)(l  * a,  ^ ) 1-1.2 

neglecting  axial  cooling.  The  constants  A,  B,  and  m were  obtained 
by  velocity  calibration  as  done  for  a single  sensor  probe;  was 
determined  by  azimuthal  yaw  calibration.  The  power  dissipated  is  pro- 
portional to  the  square  of  the  anemometer  bridge  output  minus  the  zero 
flow  voltage: 


Ei  =Mi(Q-  Qo>i  =MiUeff  (1+aiifc) 

The  voltage  can  be  linearized  directly  in  terms  of  Ueff  using  m = 
0.5.  Thus 


E*i  ■ V'/"  - MUeff  <’  + “1  U 


eff 


or  when  neglecting  higher  ordered  terms  the  instantaneous  voltage  is 
E„..  + e„4  = MUeff  [1  + r r—  + ^!r~  ] 


eff 


i U 


eff 


When  calibrated  in  a steady  flow 


Second  order  effects,  which  occur  for  angles  larger  than  20°,  are 
negligible  for  this  particular  flow. 


The  following  calibration  procedure  was  repeated  before  and  after 
each  velocity  profile  was  obtained.  The  probe  was  mounted  in  a spe- 
cially constructed  support  with  orientation  adjustment,  the  "dial 
calibrator,"  that  permitted  a ±50°  rotation  of  the  split-film  sensor 
about  a fixed  location.  The  sensor  v/as  supported  over  the  potential 
core  of  the  near  jet  produced  by  a TSI  model  1125  calibrator.  The 
dial  calibrator  was  used  for  both  the  velocity  and  azimuthal  yaw 
calibrations.  The  uncertainty  in  the  angle  was  less  than  0.2°.  The 
use  of  the  split-film  probe  requires  closely  equal  temperatures  of 
the  films,  so  that  very  fine  overheat  adjustment  of  one  film  is  man- 
datory to  match  temperatures. 

Both  films  were  exposed  to  77°  air  flow  directed  along  the  plane 
of  splits.  Cold  resistances  were  measured  and  found  to  be  about  11 
and  12  ohms.  Both  films  were  then  overheated  by  a factor  of  0.5. 

One  film  was  left  without  further  adjustment  while  the  other  film 


overheat  was  adjusted  using  a 1 ohm  precision  ten  turn  potentiometer 
The  temperature  matching  was  achieved  when  the  ratio  of  the  two  out- 
put bridge  voltages  remained  constant  (within  0.2%)  as  the  velocity 
directed  along  the  plane  of  the  splits  was  varied  from  9 fps  to  70 
fps.  Temperature  matching  was  checked  again  immediately  after  tak- 
ing data  from  flow  C in  order  to  estimate  any  possible  anemometer 
drift. 

The  bridge  outputs  were  then  each  linearized  and  checked  with 
the  known  calibrator  velocity.  A calibration  equation  of  the  form: 

7 


r 


E„ . = C.U  ,,  + D. 
ill  i eff  1 


was  determined  from  a least  squares  fit  with  less  than  4%  standard 
deviation.  The  background  noise  level  in  the  steady  laminar  calibrator 
flow  was  also  measured.  This  noise  took  the  form  of  a triangular  wave 
with  no  phase  change  between  channels.  We  suspect  the  noise  was  due 
to  vortex  shedding  off  of  the  relatively  large  diameter  sensor.  The 
probe  was  then  calibrated  for  several  azimuthal  yaw  angles  to  deter- 
mine a.j . Yaw  characteristics  were  found  to  be  independent  of  the  mag- 
nitude of  Ug^. 

The  equations  governing  the  split  film  are  a combination  of  the 
above  equations 

E(ti  + e i = ciUeff  C + “'ueff  * &i'v  + V)/Ueffl  + °1 
i « 1.  2 

’ CT  Cueff  + VJ  * D1 


or 


and 


E£2  = C2^Ueff  + 2ot2V-'  + °2 


Since  C^,  D^,  and  are  known  from  calibration,  then  and  are 
measured  and  and  V are  derived.  The  angle  of  flow,  <J>,  at  a point 
is  then  found  from 

sind>  = V/Ueff 

Similarly,  the  fluctuation  quantities  are  found  from: 

(e^)2  = C2  u2  + 4C2a-j  uv  + (20^  )2  v2 


(etl  + ei2)C  = (C1+C2)2u2+4(C1a1+C2a2)(C1+C2)uv+4(C1a1+C2a2)"v 


2, .2 


(e^  - e^2)  = (C^-C2)  u +4 ( C-jai ~C2a2 ) ( C-j -C2)uv+4 ( C-| a-j -C2a2)  v 

8 


2~2 


—2  2 2 

e 1*  ^e£.l  + e£,2^  ’ and  (e&i  " e£ 2 ^ were  measurec1’  corrected  for  the 

~Z  

above  mentioned  noise,  and  used  to  solve  simultaneously  for  u , -uv, 

~o 

and  v using  the  latter  equations.  These  results  are  discussed  in 
section  3.4  below. 

The  traversing  mechanism  used  for  the  boundary  layer  velocity 
measurements  was  mounted  on  the  supporting  frame  for  the  upper  wall 
and  provided  for  precise  positioning  of  the  probe  sensors  as  described 
by  Strickland  and  Simpson  (1973).  A cathetometer  was  used  to  accu- 
rately locate  the  probe  sensor  from  the  wall  within  an  uncertainty  of 
about  ±0.001  inches.  The  detailed  streamwise  free-stream  velocity 
distributions  were  obtained  using  a TSI  model  1210-20  hot-film  probe 
(0.002  inches  diameter  platinum-coated  quartz  rod)  attached  to  a 
probe  support  rod  inserted  from  the  tunnel  exit  or  mounted  on  the 
toy  racing  car  shown  in  figure  3 of  Simpson  and  Wallace  (1975).  The 
car  was  easily  positioned  along  the  flow  by  fishing  line.  It  could 
not  be  used  near  the  exit  since  substantial  flow  blockage  effects 
would  have  been  produced. 

The  flush  surface  mounted  hot-film  sensors  for  the  wall  "bursting" 
and  spanwise  measurements  are  described  in  detail  by  Strickland  and 
Simpson  (1973).  The  basic  sensing  part  is  a very  thin  layer  of  plat- 
inum (Engelhard  Ind.  Liquid  Bright  Platinum  #05)  fired  on  the  end  of 
a 2mm  diameter  quartz  rod.  Gold  leads  (Engelhard  Ind.  Gold  Alloy 
Paste  A- 1199)  were  fired  on  the  sides  of  the  rod  and  short  wire  leads 
were  soldered  to  the  gold.  A casing  made  from  1/4  inch  diameter  plexi- 
glas  rod  was  used  to  protect  the  sensor  from  damage  due  to  handling. 


9 


The  resulting  unit  was  mounted  in  the  wind  tunnel  wall  with  the  plati- 
num portion  flush  with  the  test  wall.  A unit  containing  two  flush 
surface  0.020  inches  diameter  platinum  sensors  was  also  fabricated 
to  permit  one  of  the  sensors  to  be  traversable.  This  unit  was  used 
in  the  surface  spanwise  structure  measurements  and  is  described  in 
more  detail  by  Simpson  and  Wallace  (1975).  All  of  these  sensors  were 
operated  at  an  overheat  ratio  of  0.03  in  the  current  experiments.  A 
higher  overheat  ratio  would  have  permanently  damaged  the  sensors. 

The  constant  temperature  frequency  response  for  each  sensor  was  deter- 
mined to  be  down  3 db  at  4 KHz  using  the  method  of  Freymuth  (1967). 

Other  electronic  equipment  included  Krohn-Hite  model  3202  and 
330B  filters,  HP  model  400E  rms  voltmeters,  an  Anadex  model  CF-600 
timer- frequency  counter,  a SAIC0R  model  41  digital  correlation  and 
probability  analyzer,  an  Applied  University  Research  four-channel  FM 
tape  recorder  (response  down  3 db  at  2 KHz),  a voltage  comparator  or 
schmitt  trigger  using  an  operational  amplifier  integrated  circuit  and 
a multiplier  using  an  Analog  Devices  AD533  JH  integrated  circuit 
trinmed  to  within  ±1%  fullscale  nonlinearity  error.  A true  integrat- 
ing voltmeter  consisting  of  a voltage-controlled  oscillator  (Tektronix 
FG501  Function  Generator)  and  a digital  counter  (Tektronix  DC503  Universal 
Counter)  was  found  to  be  superior  to  RC-type  meter  circuits  used  in 
most  voltmeters.  A HP  5451  A/71  A Fast  Fourier  Analyzer  System  was 
used  to  process  data  recorded  on  the  tape  recorder. 

For  the  intermittency  measurements  in  flow  C,  the  laser  optics, 
and  photo-multiplier  tube  of  the  SMU  laser  anemometer  (Simpson,  et  al . , 
1974)  was  used  to  scatter  light  from  the  smoke  filled  boundary  layer 
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and  collect  the  signal.  A focal  volume  of  0.32  im  diameter  and  3.56 

? 

mm  long  was  produced  from  a laser  beam  1.1  rrm  in  diameter  at  the  1/e 
intensity  locations.  Since  the  laser  anemometer  is  mounted  on  a single 
traversable  cart,  accurate  location  of  the  focal  volume  could  be 
determined. 

3.  EXPERIMENTAL  RESULTS 

3.1  Description  of  the  test  flows 

Because  of  the  strong  influence  of  acceleration  on  these  boundary 
layers,  careful  measurements  of  the  freestream  velocity  were  made 
every  one  inch  along  the  test  section  using  the  rake  and  car  probe 
mounts.  For  flow  C the  following  velocity  and  K distributions  were 
obtained: 

1/1^=  -1.525xl0'2  + (2.307x10-3)(x)  - (1  .1 33x1  O'5) (x2) 
for  96  < x < 178  inches  and 

1/U  = 4.273X10"1  - (2.466x10" 3) (x)  + (1 . 473x1 0"6 )( x2 ) 

00 

for  178  < x < 190  inches  with 

K = -4.622  x 10'6  + (4.5427  x 10'8)(x) 
for  96  < x < 185  and  K = 3.83  x 10”8  for  186  < x < 191  inches.  U is 

- - - - oo 

in  fps  and  x is  in  inches.  The  maximum  deviation  of  U data  from 
these  equations  is  less  than  0.4%  while  the  rms  deviation  is  about 
0.1%.  The  different  x ranges  for  the  1/U^  distributions  and  their 
respective  K distribution  reflect  the  fact  that  the  upstream  velocity 
equation  produces  K values  in  closer  agreement  with  data  for  1 78<x<l 85 
inches.  The  single-sample  uncertainty  (Kline  and  McClintock,  1953) 
in  K at  20:1  odd:,  is  about  2%. 
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For  flow  D,  v.ith  both  increasing  and  decreasing  K regions,  the 
following  velocity  and  K distributions  are  presented: 

1/U^  = 9.898x10" 2 + (6.568x10'5)(x) 
for  64  < x < 96  inches 

1/U  = -1 . 449x10” 1 + (5.167x10"3)(x)  - (2.675xl0"5)(x2) 

00 

for  96  < x < 132  inches 

1/U  = 3.243xlO-1  - (1.917x10'3)(x) 

oo 

for  132  < x < 144  inches 

1/U  = 8.497xlO_1  - (9.1869x10'3)(x)  + (2.51 38xlO'5) (x2) 

00 

for  144  < x < 180  inches  and 

1/U  = 2.475xl0~2  - (7.8147x10'5)(x) 

00 

for  180  < x < 190  inches. 

K = -1.32xlO-7  for  64  < x < 96  inches 

K = -1 .037x1 0-5  + (1.0737xl0“7)(x)  for  96  < x < 132  inches 
K = 3.85x10"^  for  132  < x < 144  inches 
K = 1 .8438x1 0"5  - (1.009xl0"7)(x)  for  144  < x < 179  inches 
K = 1 . 568x1 O-7  for  180  < x < 190  inches 
The  maximum  deviation  of  U data  from  these  equations  is  less  than  1 % 

oo 

while  the  rms  deviation  is  about  0.2%.  The  discontinuity  in  K at  the 
ends  of  each  region  is  of  the  order  of  10~7.  The  single-sample  uncer- 
tainty in  K is  about  4%.  The  smoothed  K distributions  for  both  flows 
C and  D are  shown  in  figure  2. 

As  was  done  for  flows  A and  B,  observations  were  made  to  assess 
the  three-dimensionality  of  flows  C and  D.  Because  of  the  care  used 
in  adjusting  the  spanwise  elements  of  the  upper  wall  parallel  to  the 
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test  wall,  the  thin  boundary  layers  in  a large  aspect  ratio  channel, 
and  the  fact  that  the  flow  was  accelerating,  minimal  three-dimensional 
effects  were  expected.  At  the  88  inches  location  there  was  less  than 
]%  spanwise  difference  in  the  momentum  thickness,  so  no  gross  three- 
dimensionality  due  to  the  upstream  flow  was  present.  The  side  wall 
boundary  layers  of  the  converging  section  for  each  flow  tended  to 
remain  at  about  a constant  thickness  due  to  the  reduction  of  the  side 
wall  surface  area  simultaneously  with  acceleration  of  the  freestream. 
Thus  convergence  or  divergence  effects  on  the  test  wall  boundary  layer 
by  the  side  wall  boundary  layers  appears  to  be  negligible. 

As  discussed  in  section  3.3  below,  the  skin  friction  was  deduced 
by  analyzing  mean  velocity  profiles  near  the  wall.  The  smoothed  "best 
estimate"  skin  friction  coefficients  and  other  required  and  experi- 
mentally deduced  quantities  were  used  to  check  the  balance  of  terms 
in  the  two-dimensional  momentum  integral  equation.  In  both  flows  C 
and  D this  equation  was  balanced  well  within  the  uncertainty  of  the 
most  uncertain  term,  de/dx,  or  about  5%  of  that  term. 


3.2  Mean  velocity  profile  measurements 

Figure  3 shows  the  mean  velocity  profile  results  for  flow  C.  A 
distinct  logarithmic  region  exists  for  each  upstream  profile,  with 
the  last  evidence  of  any  logarithmic  region  occurring  at  about  station 
160  inches.  As  shown  in  figure  4,  the  shape  factor  H is  near  a mini- 
mum at  station  165  inches.  The  wake-like  tail  of  each  profile  near 
the  freestream  is  observed  to  progressively  decrease  in  wake  strength. 
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(a)  measured  when  velocity  profile  measured. 

(b)  from  least  squares  fit  of  rake  and  cart  data. 

(c)  with  exceptions  noted,  from  faired  curve  of  data  on  figure  4 

(d)  determined  from  logarithmic  velocity  profile,  equation  (5). 


the  wake  strength  being  defined  as  the  maximum  deviation  of  the  mea- 
sured velocity  profile  from  the  extrapolated  logarithmic  region  profile. 
Downstream  of  the  location  of  the  minimum  H the  velocity  profiles  take 
on  an  increasingly  more  laminar-like  character,  although  the  stream- 
wise  fluctuation  measurements  clearly  show  that  a relatively  large 
turbulence  intensity  still  exists.  Since  the  last  5 inches  of  the 
test  flow  has  a constant  K of  about  3.83  x 10”  this  is  the  only 
region  where  constant  K asymptotic  velocity  profile  similarity  of 
f(=U/Uj  and  n(=yUOT/v)  is  possible.  In  this  region  the  velocity  pro- 
files for  n > 1000  are  similar  well  within  a deviation  of  0.01  f. 
However,  the  streamwise  fluctuation  intensity  profiles  discussed  in 
section  3.4  do  not  posses  similarity  in  this  region  and  H and  ReQ  do 
not  reach  constant  values,  as  required  for  true  asymptotic  similarity 
(Simpson  and  Wallace,  1975).  ReQ  drops  well  below  the  value  of  about 
360  (Kays  1966)  where  normal  transition  from  a laminar  to  turbulent 
boundary  layer  occurs,  so  true  laminar  behavior  could  be  expected 
downstream  had  the  test  section  been  longer  and  the  level  of  acceler- 
ation maintained. 

Flow  D also  has  a linearly  increasing  K distribution  but  only 
for  x < 132  inches.  dK/dx  is  more  than  twice  that  for  flow  C.  As 
in  flow  C,  the  initial  upstream  logarithmic  velocity  profile  pro- 
gressively decays  downstream.  Figure  5 shows  that  the  wake  strength 
also  decreases  to  zero  near  the  minimum  H location,  shown  at  about 
135  inches  in  figure  6.  The  location  of  the  minimum  H occurs  in  the 
region  where  K is  a maximum,  while  in  flow  C it  occurs  at  a much 


Figure  5.  Mean  velocity  profiles  for  flow  D.  Strea^lse  location  In  "S-re  6.  Ree>  H.  and  Cf/2  for  flow  D.  Cf/2  values  fron  table  2 

Inches  to  right  of  each  profile.  Note  displaced  ordinate  Cf  » unshifted  wall  profile;  ®»  shifted  wall  profile 

scales.  logarithmic  profile  with  K * 0.4;  □ , momentum 

Integral  equation. 
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lower  K value.  Downstream  the  profiles  assume  an  increasingly  more 
laminar-like  character  as  ReQ  drops  below  360  until  about  157.5  inches, 
even  though  K decreases  downstream  of  144  inches.  No  profile  simi- 
larity is  observed. 

Between  157.5  inches  and  about  172.6  inches  the  momentum  and 
displacement  thicknesses,  0 and  6*,  remain  about  constant,  giving  rise 
to  about  a constant  shape  factor  (H  = 6*/0 ) of  about  2.0.  As  noted 
in  section  3.7  below,  intermi ttency  measurements  downstream  of  157.5 
inches  indicate  that  considerable  turbulence  is  redeveloping.  At  the 
189.6  inches  location  a thick  logarithmic  region  500  < n < 5000  is 
observed,  so  that  this  flow  is  developing  into  a typical  low  acceler- 
ation turbulent  boundary  layer.  The  wake  strength  is  also  observed 
to  be  increasing. 

Simpson  and  Wallace  (1975)  found  that  for  each  of  flows  A and 
B the  location  of  the  H minimum  was  closely  related  to  where  the  wake 
strength  was  zero  and  where  entrainment  of  free-stream  fluid  ceased. 

For  flows  C and  D,  the  mass  flux  in  each  boundary  layer  as  reflected 
by  11^(6  - 6*)  increases  until  about  165  inches  and  135  inches,  respec- 
tively, as  shown  in  figures  7 and  8.  This  parameter  decreases  down- 
stream in  both  flows  but  again  increases  downstream  of  about  157.5 
inches  in  flow  D as  a low  acceleration  turbulent  boundary  layer  de- 
velops. These  plots  of  experimentally  deduced  values  indicate  that 
there  was  entrainment  of  non-turbulent  fluid  into  each  boundary  layer 
until  the  minimum  H locations.  As  in  the  cases  of  flows  A and  B,  the 
intermittency  measurements  discussed  below  indicate  that  after  the 
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Mass  flux  in  boundary  layer  of  flow  C 
'ftHH  line  - eauation  1211  results. 


figure  7 


x,  inches 


Figure  8 Mass  flux  in  boundary  layer  of  flow  D 
Solid  'me  - equation  (21)  results. 


v.r* 


cessation  of  entrainment  the  intermittent  turbulent/non-turbulent 
region  lies  outside  6.  Thus  there  is  no  mean  velocity  gradient  9U/3y 
and  no  wake-like  tail  in  this  region  and  no  mechanism  for  engulfment 
of  non-turbulent  fluid  by  the  intermittent  bulges  is  possible.  The 
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model  proposed  by  Simpson  and  Wallace  for  the  entrainment  process  in 
a strongly  accelerated  boundary  layer  is  discussed  in  section  4 below. 

3.3  Skin- friction  results 

The  skin-friction  coefficient  C^/2  was  primarily  determined  by 
two  methods:  the  velocity  profile  near  the  wall  and  the  logarithmic 

velocity  profile  relationship  for  unaccelerated  flow  regions.  The 
momentum  integral  equation 

+ KRee(2  + H)  (') 

was  used  some  for  flow  D when  both  terms  of  the  right  side  were  pos- 
itive and  relatively  certain. 


The  velocity  profile  near  the  wall  can  be  derived  from  the  dif- 
ferential momentum  equation,  neglecting  the  convective  and  turbulent 
transport  terms: 


Equation  (4)  was  used  with  experimental  velocity  measurements  for 
5 < y+  = n/C^/2  to  determine  C^/2.  Oka  and  Kostic  (1972)  noted  that 
hot-film  and  hot-wire  measurements  are  strongly  influenced  by  conduc- 
tion to  the  test  wall  for  y+  < 4.  For  velocity  profiles  with  loga- 
rithmic regions  it  is  known  that  equation  (4)  does  not  well  describe 
the  velocity  profile  for  y+  > 6,  so  only  points  in  the  range  5<y+<6 
were  used.  For  the  more  laminar-like  profiles  experimental  data  with 
y+  > 6 were  also  used. 

Logarithmic  velocity  profiles  in  unaccelerated  or  weakly  accel- 
erated flow  regions  are  described  by  the  relation 


where  q = G.4G  (Re  /6000)"^  (Simpson,  1970)  for  low  Reynolds  number 
0 

boundary  layers  with  Re  < 6000.  A fit  of  equation  (5)  to  experi- 
mental  data  was  made  for  the  two  upstreammost  profiles  of  flow  C and 
to  the  downs treammost  profile  of  flow  D.  In  the  latter  case  n = 0.40 
was  used  since  a large  logarithmic  region  that  normally  accompanies 
high  Reynolds  number  boundary  layers  was  observed.  Results  using 
this  latter  method  are  about  ±10%  uncertain. 

For  flow  C the  0.001  inches  diameter  cylindrical  hot-film  sen- 
sor was  located  from  the  test  wall  within  an  uncertainty  of  about 
0.001  inches.  Considering  the  ±1/2%  uncertainty  of  the  velocity,  the 
Cf/2  values  shown  in  figure  4 appear  to  be  about  ±15%  uncertain  at 
20:1  odds.  This  uncertainty  value  seems  reasonable  due  to  the  scatter 


22 


of  the  data.  The  logarithmic  velocity  profile  results  from  equation 
(5)  are  in  close  agreement  with  these  results. 

For  flow  D the  0.00015  inches  diameter  hot-wire  sensor  was  lo- 
cated from  the  test  wall  with  somewhat  more  uncertainty  than  for  flow 
C sine,?  the  wire  was  more  difficult  to  illuminate  and  to  observe  by 
the  cathetometer.  This  more  fragile  sensor  was  located  at  a safe 
distance  of  about  0.005  inches  from  the  wall  to  prevent  breaking  the 
wire.  Data  from  this  smaller  diameter  sensor  are  less  suceptible  to 
velocity  gradient  effects  so  the  near  wall  velocity  profiles  for 
y+  > 4 are  believed  to  be  of  higher  quality  than  those  of  flow  C. 

For  each  velocity  profile  point  a least  squares  fit  of  equation  (4) 
to  the  five  surrounding  data  points  was  made  to  deduce  the  slope 
9f/9n  at  that  point.  0^/2  can  be  eliminated  from  equations 
(3)  and  (4)  to  produce 


"-+/l 

9n  V 

[9fl 

9n 

1 

+ 2Kf 

K 

from  which  the  proper  sensor  distance  from  the  test  wall  could  be 
computed  from  K,  f,  and  9f/9n.  For  flow  D equation  (6)  was  used  to 
determine  the  required  change  in  the  y direction  for  a self-consis- 
tent velocity  distribution  near  the  wall.  Values  within  about  5%  of 
one  another  were  obtained  for  successive  points  in  a given  velocity 
profile.  As  given  in  Appendix  B these  shifts  in  the  y direction  are 
no  more  than  ±0.005  inches  and  therefore  are  reasonable. 

Figure  6 shows  relatively  less  scatter  in  the  wall  profile  re- 
sults than  shown  for  flow  C,  thus  an  estimated  uncertainty  of  ±10* 
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seems  reasonable.  The  logarithmic  profile  C^/2  result  at  189.6 
inches  is  in  good  agreement  with  the  C^/2  result  from  equation  (4) 
with  no  y direction  shifting.  Results  from  equation  (1)  at  161.6 
and  179.7  inches  are  also  in  good  agreement  with  the  wall  profile 
resul ts. 

Smooth  curves  of  "best  estimate"  values  for  C^/2  are  shown  on 
figures  4 and  6 with  values  given  in  Tables  1 and  2.  While  these 
curves  were  just  faired  among  the  points  shown  on  those  figures 
within  about  15%  and  10%,  respectively,  they  are  probably  represen- 
tative of  the  variation  of  C^/2  along  each  flow.  From  the  laminar 
sink  flow  (Schl ichting,  1968),  i.e.  a constant  K flow,  one  can  obtain 
the  expression 


so  for  the  maximum  K of  about  3.83  x 10  ^ achieved  in  these  two  flows, 

- "5 

C^/2  = 2.26  x 10  . Near  the  maximum  K for  each  of  flows  C and  D 

the  experimental  C^/2  values  are  above  this  value  as  they  should  be 
since  considerable  turbulent  momentum  transport  is  still  present. 

While  one  cannot  exactly  compare  results  from  different  K distribu- 
tions, it  is  interesting  that  the  experimental  results  shown  in  fig- 
ure 6 for  flow  D are  always  greater  than  or  about  equal  to  those 
given  by  equation  (7)  with  the  local  K value.  In  the  region  near 
165  inches  where  the  turbulence  intermi ttency  is  near  zero  at  the 
wall,  as  discussed  in  section  3.7  below,  the  skin  friction  values 
are  about  equal.  Downstream  of  this  region  the  experimental  C^/2 
values  increase  rapidly  as  turbulence  redevelops. 
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3.4  Reynolds  stresses  distributions  

FT 

Figures  9 and  10  show  intensity  v u /U^  profiles  for  flows  C 
and  D from  the  fluctuation  portion  of  the  linearized  hot-film  and 
hot-wire  signals.  In  flow  C the  maximum  intensity  remains  at  about 
0.1  and  is  located  in  the  range  of  300  < n < 400  up  until  K>3.6xl0~^ 
at  about  181  inches,  after  which  the  normalized  intensity  decays  and 
the  maximum  intensity  location  moves  to  a greater  n value.  The 
hump  in  the  outer  region  of  the  upstream  intensity  profiles  vanishes 
just  downstream  of  the  cessation  of  entrainment  at  about  165  inches. 
Between  165  inches  and  181  inches  the  acceleration  of  the  wall  region 
flow  alone  is  apparently  sufficient  to  maintain  the  same  level  of 
turbulence  production  as  upstream.  In  other  words  the  entrainment 
of  high  momentum  free-stream  fluid  performs  no  role. 

Flow  D exhibits  the  same  behavior  but  the  outer  region  hump  does 
not  completely  vanish  until  about  143.8  inches.  The  maximum  inten- 
sity is  located  in  the  range  of  200  < n < 400.  Beginning  at  about 
153.6  inches,  where  C^/2  begins  to  drop  rapidly,  the  value  and  n 
location  of  the  maximum  intensity  decreases.  After  retransition  to 
a low  acceleration  turbulent  boundary  layer  begins, the  maximum  inten- 
sity remains  at  about  0.08  but  the  n location  of  the  maximum  increases. 
At  189.6  inches  the  intensity  profile  of  a normal  low  acceleration 
turbulent  boundary  layer  is  emerging. 

~2  ~2  

The  split-film  probe  was  used  to  determine  u , v , and  -uv  for 
flow  C.  No  split-film  probe  measurements  were  attempted  for  flow  D. 

The  mean  velocity  profiles  showed  very  good  agreement  with  those 
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g ure  10.  Stream!  $e  intensity  profiles  for  flou  0.  Note  displaced 

scales.  Station  position  in  Inches  at  right  of  each  profile, 
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measured  with  a single  film  boundary  layer  probe.  The  deviation  was 
within  the  uncertainty  of  the  calibration  except  at  locations  171.3, 
178.8,  and  187.6  inches.  At  these  locations  the  mean  velocities 

3 

below  n = 10  were  as  much  as  47%  high  nearest  the  wall  yet  were  less 
than  2%  high  at  n=10^  and  were  less  than  0.5%  high  at  n = 1300  and 
beyond.  At  locations  101.7,  137.8,  149.6  and  162.4  inches,  the  mean 
velocity  deviated  from  the  single  film  probe  by  less  than  3%  from 
the  wall  to  the  freestream.  The  large  deviation  near  the  wall  at 
higher  freestream  velocities  may  be  due  to  the  significantly  greater 
velocity  gradient  effects  which  tend  to  cool  the  top  film  more  than 
the  bottom  film. 

The  V component  of  velocity,  when  normalized  against  U ^ allowed 
comparison  of  the  angle  of  flow  <f>  with  respect  to  the  test  wall  since 
sin  4>  = V/U^p.  After  accounting  for  the  angle  of  the  plane  of  the 
splits  with  the  test  wall,  the  angle  of  flow  measured  by  the  split- 
film  probe  at  a point  in  the  freestream  was  at  worst  within  1.4°  of 
the  flow  angle  at  that  point  derived  from  the  displacement  thickness 
gradient  and  freestream  velocity  gradient.  This  difference  is  beyond 
0.6°  uncertainty  in  computing  <}>  from  the  calibration  equations  plus 
the  0.2°  uncertainty  in  <p  due  to  approximation  of  the  displacement 
thickness  gradient.  However,  for  6 of  the  8 profiles  taken,  the 
angles  measured  agree  within  that  0.6°  uncertainty  in 

The  u fluctuations,  as  seen  in  figure  11,  tend  to  show  very 
poor  agreement  with  those  measured  with  the  single  film  probe  for 

3 

n < 10  . At  x = 137.8  and  149.6  the  discrepancy  is  apparent  to 
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Figure  11.  Streamwise  fluctuation  distributions  from  the  snlit-filn 

sensor  for  flow  C.  Profile  location  in  inches  at  right  of 
each  profile.  Dashed  lines  - normal  hot-film  results  from 
figure  9 for  nearby  location.  Note  displaced  ordinate 
scales. 
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q = 10  . In  all  other  cases  the  agreement  is  good  from  n = 10  to 

the  freestream  except  at  x = 158.3  where  agreement  is  poor  until  n 
reaches  5000.  The  kinematic  Reynolds  shearing  stress  -uv  also  shares 

3 

the  n = 10  cutoff  above  which  the  data  begin  to  seem  reasonable. 

There  is  a rapid  decrease  in  -uv  as  the  flow  accelerates  as  shown 
in  figure  12.  Also  shown  on  this  figure  is  the  shearing  stress  pro- 
file for  190.1  inches  generated  from  the  momentum  equation  and  U/U^ 
vs.  n similarity.  The  outer  region  mean  velocity  profiles  at  185.8 

and  190.1  inches  are  similar  well  within  0.01  in  U/U  and  K is  con- 

00 

stant.  Thus  the  similarity  assumptions  are  met  for  an  asymptotic 
flow  and  the  equation 

^ = K I (1  - f2)dn  - (8) 

lr  Jq  drl 

describes  the  Reynolds  shearing  stress  (Simpson  and  Wallace,  1975). 

The  agreement  must  be  considered  good  for  q > 2600.  Thus  the  results 

from  the  split-film  for  -uv  are  credible  at  least  in  locations  where 
~2  ~~2 

u is  credible.  The  v quantities  measured  by  the  split-film  are 
high  in  the  viscous  sublayer  due  to  the  high  velocity  gradient  there. 
The  data  in  the  outer  region,  q > 3000,  are  more  believable. 

Many  problems  were  encountered  in  using  the  split-film  sensor. 

The  calibration  process  was  time-consuming  and  frequent  checks  were 
required  to  overcome  anemometer  drift.  Anemometer  drift  at  times 
caused  the  films  to  operate  with  different  temperatures,  thus  produc- 
ing thermal  interference  between  films  as  well  as  introducing  a 
slight  velocity  dependence  on  the  yaw  characteristics.  The  resistance 


29 


of  the  probe  support  structure  sometimes  changed  when  the  probe  was 
unplugged  and  then  replaced  in  its  support.  Vortex  shedding  from 
the  sensors  produced  high  noise  levels  at  increasing  velocities. 

This  noise  introduced  greater  uncertainty  in  the  fluctuation  quanti- 
ties. Velocity  gradients  near  the  wall  produced  unusually  high  V 
values  at  those  points  where  V should  be  lowest.  Several  schemes 
were  tried  to  correct  for  the  influence  of  the  high  velocity  grad- 
ient. Unfortunately,  those  corrections  did  not  produce  consistent 
results. 

This  experience  with  the  split-film  probe  seems  to  be  reason- 

I 

ably  consistent  with  that  of  Sandborn  (1976).  He  found  v u /U^ 
to  be  closely  the  same  as  that  measured  by  a single  hot  wire  paral- 
lel to  the  test  wall,  -uv  values  from  split- film  data  appeared  to 
be  close  to  the  estimated  true  shear  stress  profile  away  from  the 
wall  where  velocity  gradients  are  sufficiently  small.  On  the  other 
hand.  Young  (1976)  found  that  the  split-film  probe  consistently  pro- 
duced a low/7  by  about  22%  and  a low  -uv  by  about  26%  in  a two- 
dimensional  channel  flow.  Thus,  clearly  the  split-film  results  should 
be  critically  compared  with  results  from  other  sensors  for  consistency. 


3.5  Spectra  measurements  and  dissipation  estimates 

The  spectrum  function  F(n)  of  u was  obtained  throughout  each 
boundary  layer  at  several  streamwise  locations,  where 


F(n)dn  = u^ 


(9) 
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and  n is  the  frequency  in  Hz.  Simpson  and  Wallace  (1975)  discovered 
that  a portion  of  F ( n ) varied  like  n ^ for  locations  in  the  region 
n > 1000  after  the  cessation  of  entrainment  in  flow  B.  They  found 
that  Tchen's  high  mean  vorticity  spectral  model  (Hinze,  1975)  relat- 
ing F(n)  to  the  turbulence  dissipation  e by 

~~2 

£*=  |a[nF(n)]  % §£  (10) 

U U n 

OO  CO 

seemed  to  hold. 

Their  data  indicate  that  a is  about  0.77  for  flow  B.  Data  from 
flow  B were  reanalyzed  using  equation  (10)  above.  Equation  (21)  of 
Simpson  and  Wallace  should  have  contained  u . The  spectral  data  of 
Laufer  (1954)  in  the  logarithmic  region  of  a pipe  flow  indicate  that 
a ~ 0.88  while  Klebanoff's  (1955)  flat  plate  turbulent  boundary  layer 
produced  flat  nF(n)  spectra  with  a ~ 0.80.  In  both  of  these  latter 
cases,  e was  evaluated  by  the  equilibrium  relationship 


e = -uv 


that  is  applicable  for  the  logarithmic  velocity  profile  region  (Rotta, 
1962).  Here  we  wish  to  examine  the  spectra  obtained  for  flows  C and 
D and  the  range  of  applicability  of  equation  (10)  using  a = 0.8. 
Representative  first  moments  of  the  spectra  nF(n)  for  flows  C 


and  D are  presented  in  figures  14  through  23.  The  hot-film  and  hot- 
wire sensor  response  was  down  5%  at  15  KHz,  and  25  KHz,  respectively. 


The  data  were  recorded  on  tape  and  processed  on  the  fast  Fourier 


analyzer  with  a bandwidth  of  2 Hz  and  50  seconds  **ecord  times.  The 
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2 KHz  upper  frequency  limit  on  the  analyzer  was  not  a serious  limi- 
tation except  for  the  data  at  the  downstream  stations  in  flow  D. 

For  several  of  the  spectral  distributions,  a small  range  of  the 
-5/3  law  of  the  inertial  subrange  existed  as  shown  in  figure  14,  so 
it  was  used  to  estimate  the  dissipation  rate  e: 


U2  F(n)  - KE2'3 


f 2ir 

U 


^-2/3 


■5/3 


02) 


Here  K is  a constant  taken  to  be  0.49  (Corrsin,  1 9G4 ; Bradshaw, 

1967a)  and  U is  the  local  mean  velocity.  Bradshaw  (1967b)  suggested 

-j  1/2 

that  the  turbulence  Reynolds  number  Re^  = (u  ) A/v,  where 
2 ~2 

A = 15v  u /e,  must  be  greater  than  100  for  an  inertial  subrange  to 
exist.  For  those  spectra  on  which  a region  with  a -5/3  slope  was 
found,  35  < Re  < 80.  Thus  there  is  some  question  of  the  validity 

A 

of  equation  (12)  even  for  these  cases. 

Figures  24  and  25  show  the  dissipation  rate  results  for  flows 
C and  D.  For  comparison  the  dissipation  rate  for  a turbulent  bound- 
ary layer  which  possesses  a logarithmic  "law  of  the  wall"  mean  velo- 
city profile  is  given  by 
f C, 


ev  _ 
U4 


2 

2 


3/2 


1 


0.4ri 


(13) 


This  equation  is  derived  from  equation  (11).  This  result  would  also 
be  expected  to  be  valid  for  moderately  accelerated  turbulent  bound- 
ary layers  since  a logarithmic  mean  velocity  profile  would  still 

4 

exist.  At  the  upstream  stations  the  dimensionless  dissipation  ev/U^ 
is  the  same  order  of  magnitude  for  each  flow  from  equations  (10), 
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(12),  and  (13).  The  results  from  equation  (12)  for  each  flow  near 

the  100  inches  location  are  closely  the  same.  At  137.2  inches  in 

flow  C the  results  from  equations  (10)  and  (12)  are  reasonably  close 

for  a given  n when  a = 0.8  is  used.  Since  the  nF(n)  spectra  at 

downstream  stations  are  flat  for  n > 800,  the  results  from  equation 

(10)  are  the  best  available  estimates  of  the  dissipation  rate. 

Spectral  data  from  locations  170.8,  178.3,  182.3,  and  185.8 

inches  of  flow  C are  very  similar.  In  the  sublayer,  n < 800,  no  flat 

region  of  a given  nF(n)  distribution  is  observed.  The  contribution 

from  the  20  Hz  free-stream  unsteadiness  becomes  an  increasingly 

“2 

greater  portion  of  the  total  u in  the  downstream  direction.  Figure 

9 shows  that  the  total  contribution  to  u"  at  a given  n decays  down- 

stream  of  about  178  inches.  However,  since  (u  ) (nF(n) ) is  the  same 

whether  the  unsteadiness  is  accounted  for  or  not,  the  results  in 

figure  24  should  be  unaffected.  For  a given  n in  the  outer  region, 

4 

ev/U^  decreases  two  orders  of  magnitude  from  137.2  inches  to  170.8 
inches.  Thereafter  the  decrease  is  much  slower.  In  the  inner  region, 

4 

n < 1000,  cv/U^  remains  about  the  same  order  of  magnitude  all  along 
the  flow. 

Figure  19  shows  spectral  distributions  at  135.6  inches  for  flow 
0.  This  location  is  the  last  downstream  station  that  was  recorded 
with  spectra  shapes  similar  to  those  for  flow  C.  Between  this  loca- 
tion and  161.6  inches,  retransition  to  an  entraining  turbulent  bound- 
ary layer  begins,  as  mentioned  in  section  3.7  below.  Figure  20  shows 
that  F(n)  varies  as  n 


in  the  region  of  r,  < 800  at  161.6  inches, 


where  the  intermittency  is  about  0.01.  In  the  outer  region  at  this 
location  flat  nF(n)  spectra  are  found.  Note  the  relative  importance 
of  the  20  Hz  unsteadiness.  Downstream  where  the  intermittency  is 
progressively  increasing,  the  spectral  distributions  are  flatter. 

At  179.7  and  189.7  inches  the  peaks  shift  progressively  toward  1 KHz 
as  the  low  accelerated  turbulent  boundary  layer  redevelops.  Notice 
that  20  Hz  unsteadiness  does  not  dominate  the  near  wall  spectra. 

The  dissipation  results  for  flow  D from  equations  (1C)  and  (12) 
1 are  in  reasonably  good  agreement  at  123.1  inches  and  in  fair  agree- 

ment at  132.2  and  140.7  inches.  Between  132.2  and  161.1  inches  the 
results  from  equation  (10)  are  considered  the  best  available  esti- 
mates of  the  dissipation  rate.  For  a given  n in  the  outer  region, 

4 

ev/U^  decreases  two  orders  of  magnitude  from  123.1  to  161.1  inches, 

4 

as  was  observed  for  flow  C.  As  for  flow  C,  ev/U  remains  about  the 

00 

same  order  of  magnitude  in  the  inner  region  along  this  length  of 
flow.  Certainly  the  spectral  behavior  downstream  of  retransition 
needs  further  research. 

3.6  Wall  bursting  frequencies  and  spanwise  structure 

The  more  or  less  periodic  1 i Ft  off  or  bursting  process  in  the 
viscous  sublayer  has  become  accepted  as  the  sequence  of  events  tnat 
create  the  Reynolds  shearing  stress  near /the  wal  l (Wallace  et  al . , 
1972;  Nychas  et  al . , 1973;  Willmarth,  1^975).  On  a short-time  basis 

i 

the  velocity  fluctuations  in  the  sublayer  produced  by  the  bursting 
process  vary  across  the  sublayer  in  the  spanwise  direction.  The 
data  of  Gupta  et  al . (1971)  indicate  that  the  short-record-time 

i • 
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normalized  cross-correlation  of  fluctuations  sRuu(z,  Ts)  could  be 
represented  by  a spanwise  periodic  function.  Their  results  indicate 
that  the  short-record-time  T$  over  which  the  cross-correlation  is 
averaged  must  be  less  than  about  20  bursting  periods  in  order  for 
the  periodic  spanwise  structure  to  be  detectable.  Simpson  and 
Wallace  (1975)  measured  the  average  frequency  of  bursting  and  the 
long-time  spanwise  cross-correlation  of  fluctuations  for  two  sink- 
flow  type  strongly  accelerated  turbulent  boundary  layers,  flows  A 
and  B.  Here  the  results  from  flows  C and  D are  presented. 

Strickland  and  Simpson  (1973,  1975)  assumed  that  the  short-time 
autocorrelation  time  scale  from  a flush  wall  hot-film  sensor  and  the 
bursting  period  were  equal.  Histograms  of  the  frequency  character- 
ized by  the  time  to  the  first  peak  on  these  autocorrelations  were 
constructed  and  the  corresponding  frequency  of  the  peak  of  each 
histogram  was  taken  as  the  characteristic  frequency.  The  histograms 
appeared  to  have  a log-normal  probability  distribution,  so  this  oeak 
frequency  was  also  the  median  frequency.  They  also  proved  a one  to 
one  correspondence  between  this  characteristic  frequency  and  the 
peak  of  the  first  moment  of  the  wall  shearing  stress  spectra  nF(n). 
Also,  examination  of  spectral  data  from  zero  pressure  gradient  bound- 
ary layers  produced  bursting  frequencies  in  agreement  with  those 
previously  reported.  More  recently,  Simpson  (1976)  used  a pattern 
recognition  algorithm  to  determine  the  average  bursting  frequency 

i 

from  flush  surface  hot-film  signals  produced  in  the  Max-Planck- 
Institut  fUr  Strijmungsforschung  oil  channel.  This  average  bursting 
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frequency  closely  agreed  with  the  peak  of  the  nF(n)  spectral  distri- 
bution of  these  signals. 

As  was  done  for  flows  A and  B,  the  peak  of  each  nF(n)  curve  for 
flows  C and  D was  deduced  as  the  bursting  frequency  for  a given  wall 
spectral  distribution.  The  results  are  given  in  table  3.  As  i t is 
clear  from  figures  26  and  27,  it  is  difficult  to  select  a precise 
single  frequency  at  which  nF(n)  is  a maximum  for  a given  distribu- 
tion. Consequently,  a range  over  which  the  peak  frequency  definitely 
occurs  is  also  presented  in  table  3. 

Two  types  of  non-dimensional izing  parameters  have  been  used  in 
previous  bursting  frequency  studies,  inner  variables  and  outer  vari- 
ables. Inner  variables  are  the  shear  velocity  for  the  velocity 
scale  and  v/U^.  for  the  length  scale.  This  scaling  does  not  corre- 
late other  available  bursting  data.  The  bulk  of  available  bursting 
frequency  results  for  turbulent  boundary  layers  correlate  using  the 
outer  flow  velocity  and  length  scales,  Uro  and  6.  For  zero  pressure 
gradient  boundary  layers  IJ^/Sn^  is  about  5 (Rao  et  al . , 1971)  and 
varies  between  11.7  and  8,35  for  the  separating  turbulent  boundary 
layer  of  Simpson  et  al . (1974).  It  should  be  noted  that  this  large 
eddy  outer  variables  scaling  continued  even  after  the  boundary  layer 
separa ted. 

For  flow  A presented  by  Simpson  and  Wallace  (1975),  U^/ 6nb 
started  at  about  5 at  89.2  inches,  decreased  to  about  3.5  at  about 
139  inches  or  the  location  of  minimum  H,  and  then  monotonically 
increased  to  a value  above  50  at  the  last  measurement  station.  The 
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bursting  frequency  nb  apparently  reached  about  a constant  value  after 

the  asymptotic  condition  was  achieved,  so  an  increasing  Ura/5nb 

reflects  the  fact  that  U increases  and  6 decreases  faster  than  n. 

o°  b 

increases.  For  flow  B U^/ 5 n b remained  nearly  constant  at  about  3(±1) 
until  the  location  of  a minimum  in  H,  downstream  of  which  it  also 
increased  monotonically  to  a value  of  about  90  at  the  last  measure- 
ment station.  The  bursting  frequency  n^  was  virtually  constant  ail 
along  flow  B.  The  outer  variables  parameter  UOT/6*n^  was  no  better 
a correlation  parameter  for  these  data. 

The  results  for  flow  C indicate  that  the  bursting  frequency 
increased  up  to  about  165.3  inches  or  where  H was  about  a minimum 
and  where  entrainment  ceased.  The  plateau  on  the  nF(n)  spectral 
distributions  existed  over  an  increasingly  wider  frequency  range  at 
downstream  locations.  The  low  frequency  end  of  the  plateau  remained 
at  about  25  Hz.  The  results  from  these  three  flows  indicate  that 
they  are  non-equilibrium  flows  at  least  up  until  the  cessation  of 
entrainment.  Thus  no  equilibrium  non-dimensional  parameter  such  as 
uySn^  should  have  a constant  value  upstream  of  the  minimum  in  H. 
Downstream  of  the  minimum  H location  there  was  no  entrainment,  so 
the  large  eddy  structure  and  the  intermittent  turbulent/non-turbulent 
outer  region  flow  were  not  creating  new  bulges.  Thus  the  bursting 
frequency  and  the  frequency  of  passage  of  the  intermittent  bulges 
(see  section  3.7  below)  were  constant  in  this  downstream  region. 

For  flow  D the  bursting  frequency  increased  until  near  the  loca- 
tion where  H was  a minimum,  after  which  it  remained  constant  until 
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Flow  C 

Station 

(inches) 

Peak 

frequency, 

v Hz 

Range  of 
peak,  Hz 

90.2 

105.0 

110.4 

117.3 

122.0 

129.9 

131 .6 

135.5 

138.6 

139.7 

147.6 

156.3 

158.3 

165.3 

171.8 

175.2 

177.4 
1 82.0 

184.2 

186.2 

Table  3. 


Station 


13-18 

10-15 


13- 20 
11-20 
11-20 

14- 22 
10-20 
10-20 

15- 30 
15-20 
15-25 


15-35 

20-30 

20-30 

20-45 

20-100 


Station 
(inches ) 


105.0 

110.4 

122.0 
129.9 

133.5 

135.5 

139.7 

147.6 

158.3 

171.8 

175.2 

177.4 
182.0 

184.2 

186.2 

189.6 


Flow  D 


Peak 

frequency, 
n.  , Hz 


Range  of 
peak,  Hz 


250- ?50 
150-250 
400-700 
400-700 
400-700 


Peak  frequencies  of  nF(n)  spectral  distributions  from  flush 
surface  hot-film  sensors. 


2t  . 
mm-max 


, + 

XZ2 

L^.mrn 

LZ 

156.3 

165.3 
182.0 


10.16 

9.40 

8.38 


99  4.88 

151  4.46 

142  3.23 

119  1.42 


14.48  179 
9.40  160 
8.38  274 
7.87  332 


2.41 
3.30 
21  .6 


12.50 


33.3 

44.3 

36.4 
152.5 

1007 


Table  4.  Experimental  results  from  the  wall  spanwise  spatial  cross 
correlatives.  First  five  entries,  flow  C;  second  five, 
flow  D. 


the  beginning  of  retransition  to  a low  accelerated  turbulent  bound- 
ary layer.  As  discussed  in  section  3.7  below,  high  velocity  pulses 
are  first  detected  at  about  165.3  inches.  Figure  27  shows  that  the 
wall  sensor  spectra  are  not  appreciably  affected  until  about  175.2 
inches. 

Figure  28(c)  shows  a typical  signal  trace  for  this  location. 
Between  175.2  and  182.0  inches,  a distinct  local  peak  is  observed 
in  the  spectra  at  frequencies  an  order  of  magnitude  higher  than 
observed  upstream.  This  indicates  that  the  high  frequency  oscilla- 
tions that  are  contained  in  the  higher  velocity  pulses  are  increas- 
ingly important.  Downstream  of  184.2  inches  the  laminar-like  low 
velocity  regions  such  as  seen  in  figure  31  for  182.0  inches  were  no 
longer  observed  on  signal  traces.  Thus  the  flow  at  the  wall  was 
always  turbulent.  Peak  frequencies  are  observed  in  figure  27  to 
occur  at  an  even  higher  frequency  of  about  500  Hz. 

Several  researchers  have  proposed  that  during  the  bursting  pro- 
cess a hairpin  or  horseshoe  vortex  with  trailing  legs  is  formed 
(Willmarth,  1975).  With  this  model,  the  counter- rotatinq  trailing 
legs  would  produce  fluctuations  in  the  sublayer  that  are  periodic 
in  the  spanwise  or  z direction  on  a short-time  basis.  Simpson  (1976) 
noted  that  this  model  suggests  motions  too  coherent  to  be  consistent 
with  experimental  observations. 

He  alternatively  proposed  that  "fingers"  of  high  velocity  fluid 
come  from  the  outer  region  at  a velocity  much  higher  than  sublayer 
mean  velocities  toward  the  wall  and  displace  the  low  velocity  fluid 
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top  trace-movable  sensor,  500  Hz  low  pass  filtered,  lower  trace-fixed  sensor,  1 kH 
low  pass  filtered;  10  msec/div.  Middle  photo:  171.5  inches,  20  msec/div.  Bottom 
photo:  175.0  inches,  5 msec/div. 


t 
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gure  29.  Flush  surface  hot-film  signals  from  two-sensor  unit  at  156.3  inches  for  flow  D.  Top 

trace  of  each  photo  - moveable  sensor;  bottom  trace  - fixed  sensor.  From  top  to  bottom 
Z = 0 mm,  2.29  mm,  and  8.38  mm.  Abscissa  - 10  msec/div.  Low  pass  filter:  500  Hz 

moveable  sensor;  2 kHz  - fixed  sensor.  G denotes  mean  signal  level. 


t -> 

Figure  30.  Flush  surface  hot-filn  signals  from  two-sensor  unit  at 


165.3  inches  for  flow  D.  Top  trace  of  each  photo  - 
moveable  sensor;  bottom  trace  - fixed  sensor.  Abscissa  - 
10  msec/div.  low  pass  filtered  at  2 kHz.  G denotes  mean 
signal  1 evel . 
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trace  of  each  photo  - fixed  sensor;  bottom  trace  - moveable  sensor.  From  top  to  bottom 
AZ  = 0 mm,  8.13  mm,  and  19.56  mm.  Abscissa  - 5 msec/div.  Low  pass  filter:  1 kHz  - 

moveable  sensor;  5 kHz  - fixed  sensor.  G denotes  mean  signal  level. 


to  each  side.  These  fingers  of  high  velocity  fluid  were  found  for 
the  MPI  oil  channel  flow  to  be  no  less  than  Az+  = 11  nor  greater 
than  about  az+  = 44  in  size  and  are  spaced  spanwise  across  the  flow 
at  average  distances  of  about  X^IOO.  There  are  curved  fronts  on 
these  higher  velocity  fluid  fingers  so  that  the  low  velocity  fluid 
just  downstream  at  that  instant  must  move  aside  to  satisfy  continu- 
ity requirements.  Thus  low  velocity  fluid  trapped  between  two  high 
velocity  fluid  fingers  must  move  outward  from  the  wall. 

After  the  low  velocity  fluid  has  been  "ejected"  away  from  the 
wall,  the  two  adjacent  higher  velocity  fluid  fingers  coalesce.  Since 
this  fluid  now  occupies  the  wall  region  it  progressively  gives  up 
its  momentum  until  new  high  velocity  fingers  from  the  outer  region 
force  this  fluid  away  from  the  wall  and  the  process  is  repeated  down- 
stream. The  spanwise  locations  of  these  high  velocity  fingers  vary 
randomly  for  successive  burst  occurrences. 

An  interesting  point  is  that  these  high  velocity  fingers  must 
be  formed  so  that  the  low  velocity  fluid  can  be  ejected  between  them. 
Otherwise  a blanket  of  higher  velocity  fluid  would  trap  the  low  velo- 
city fluid  beneath.  Stability  considerations  then  require  regions 
of  high  velocity  fluid  separated  by  low  velocity  ejections.  The 
bulk  of  available  experimental  results  for  low  pressure  gradient  flows 
show  that  the  most  preferred  spacing  of  these  higher  velocity  regions 
is  X*  = 100.  It  appears  unlikely  that  streamwise  vortices  that  ro- 
tate more  than  one  revolution  are  produced  by  the  ejection  process. 

Some  characteristics  of  the  spanwise  sublayer  spatial  structure 
for  flows  C and  D were  determined  using  the  two  sensor  wall  unit 
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described  by  Simpson  and  Wallace  (1975).  The  unit  was  located  span- 
wise  across  the  tunnel  with  the  direction  of  travel  of  the  slider 
plate  perpendicular  to  the  streamwise  flow  direction.  The  signals 
from  the  two  sensors  were  time-delay  cross-correlated  using  the 
SAICOR  model  41  correlator,  with  the  record  time  for  a given  z-spac- 
ing  time-delayed  correlation  being  65.5  seconds.  For  zero  time 

A 

delay,  normalized  spanwise  cross-correlations  sR^(z)  such  as  shown 
in  figures  32  and  34  result.  Three  characteristic  lengths  were 
deduced,  the  distance  to  the  first  correlation  maximum  z . the 

ITlaX 

distance  to  the  first  correlation  minimum  z . , and  the  integral 

mi  n 

length  scale 


Lz  = 2 


z . 
1 


sR(z)dz 


(14) 


The  limit  of  integration  z^  was  taken  as  the  largest  available  loca- 
tion for  which  data  were  available. 

Figure  32  and  table  4 indicate  that  for  flow  C the  integral 
length  scale  progressively  decreased  until  165.3  inches  and  slightly 

increased  thereafter.  The  parameters  z . and  z are  relatively 

min  max 

easy  to  interpret,  within  1 0%  for  the  former  and  20"  for  the  latter. 
The  results  of  Gupta  et  al . (1971 ) and  Simpson  (1  975)  indicate  that 
zmin  a somewhat  more  reliable  spatial  parameter  than  zmax-  These 
two  parameters  are  normalized  by  the  wall  length  scale  v/U^  to  pro- 
duce the  quantities 


, U z , 2U  z . 

,+  _ t max  __j  ,+  _ t min 

A,-,  = and  \ 0 — ■ — 

zi  v z d v 


(15) 
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which  have  the  value  of  about  100  for  zero  pressure  gradient  bound- 
ary layers.  Only  at  182.0  inches  does  A^  increase  appreciably  from 

A 

this  value.  At  the  two  downstream  locations,  sR^r(z)  behaves  more 
like  the  low  Reynolds  number  cross-correlation  computed  by  Simpson 
(1976).  He  used  the  Gupta  el  al . spanwise  spacing  probability  dis- 

A 

tribution  to  compute  sRtt(z)  when  no  large  eddy  outer  region  struc- 
tures were  present,  i.e.  L /z  -n  -*■  0. 

Figure  34  and  table  4 indicate  the  same  type  of  behavior  for 
flow  D up  to  156.3  inches  with  Lz/znn-n  ~ 1/2.  Figures  28(a)  and  29 
shows  typical  simultaneous  signal  traces  for  the  two  sensors.  For 
the  two  downstream  stations,  the  signals  contained  high  velocity 
pulses  intermittently  as  shown  in  figures  30  and  31  and  as  discussed 
in  section  3.7  below.  The  length  scale  ratio  Lz/Zm^n  was  about  1.1 
at  165.3  inches  and  0.579  at  182  inches.  The  parameter  A^  is  sig- 
nificantly larger  than  100  at  these  two  stations. 

Several  different  models  were  examined  in  an  attempt  to  corre- 
late these  Az2  results.  As  attempted  for  flows  A and  B,  A^  was 
compared  with  the  data  of  Schraub  and  Kl i ne( 1 965)  when  presented 
versus  K.  Just  as  Schraub  and  Kline  found,  A^  is  of  the  order  of 
300  for  K = 3 x 10  However,  for  lower  K values  A^  varied  rather 
randomly  between  100  to  164.  In  flow  D A^  remained  only  at  about 
158  at  156.3  inches  even  after  the  flow  had  been  subjected  to  a K 
of  3.85  x 10  6 and  K was  decreasing.  After  retransition  to  a low 
accelerated  turbulent  boundary  layer  had  begun  A ^ was  very  large 
even  though  K was  of  the  order  of  10"®.  It  does  not  aopear  that  A^ 
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is  solely  a function  of  K.  After  retransition  began,  the  flow  char- 
acter with  high  velocity  pulses  is  entirely  different  than  that 
found  upstream. 

Simpson  et  al . (1974)  found  that  for  adverse  pressure  gradient 

turbulent  boundary  layers,  the  maximum  turbulent  shearing  stress 

1 / 2 

should  be  used  in  the  velocity  scale  in  order  to  produce  (-uv  ) ' 

ITIa  X 

z /v  values  of  about  100.  This  model  was  not  successful  for  flows 

Ilia  X 

A and  B (Simpson  and  Wallace).  Unfortunately  the  split-film  results 

for  -uv  presented  in  section  3.4  are  not  reliable  in  the  wall  region 

— 2 

where  it  is  a maximum.  Since  -uv/U^  decreases  along  a strongly 

I /O 

accelerated  flow,  (-uvm  ) z /v  would  also  decrease.  This  para- 

ITlaX  niaX 

meter  would  fall  well  below  100  since  (-uvm  )/tw  w°uld  be  below  1/4 

2 -1 12 

at  say,  156.3  inches.  Another  correlation  parameter,  (KU^/vn^) (C^/2) 

used  by  Simpson  et  al . was  also  not  successful  in  correlating  these 
data. 

It  appears  that  the  A*  spacing  lags  behind  local  conditions. 

In  other  words,  the  spatial  parameters  z and  z . remain  large 

max  mi n 

due  to  upstream  conditions  in  these  non-equilibrium  flows.  For 
example,  at  182.0  inches  in  flow  C,  U^/v  is  more  than  twice  that  at 
165.3  inches  but  z . and  z are  only  slightly  different.  In 
retrospect,  there  is  no  reason  to  expect  flows  C and  D to  possess 
equilibrium  characteristics. 

Figures  33  and  35  show  the  local  maxima  and  minima  from  the 
time-delayed  cross-correlations  for  a given  z soacing  which  were 
clearly  distinguishable  from  noise.  The  wavelike  nature  of  the 
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viscous  sublayer  behavior  is  illustrated  by  the  fact  that  the  maxi- 
mum correlation  at  a given  z spacing  is  time-delayed.  The  bursting 
frequency  results  given  in  table  3 above  are  somewhat  crudely  sup- 
ported by  the  frequency  corresponding  to  twice  the  time  delay  between 
the  first  minima  and  the  first  maxima,  as  given  in  table  4.  In  other 
words,  if  there  is  a repetitive  wave-like  nature  of  the  sublayer  flow 
in  both  z and  t,  the  time-averaged  characteristic  bursting  frequency 
should  be  approximately  the  same  at  all  z sensor  spacings. 

For  flow  D figures  29(b)  and  30(a)  show  that  for  the  z spacings 
of  sensors  at  which  the  spanwise  cross-correlations  were  a minimum, 
the  two  signals  were  strongly  out  of  phase.  At  time  A in  the  latter 
oscilloscope  photograph  a short  period  of  high  velocity  fluid  passed 
over  the  moveable  sensor.  At  the  same  time  lower  velocity  fluid 
passed  over  the  fixed  sensor.  A short  time  later,  of  the  order  of  3 msec, 
high  velocity  fluid  passed  over  the  fixed  sensor  while  low  velocity 
fluid  passed  over  the  moveable  sensor.  About  15  msec  later  the  two 
traces  again  appear  to  be  strongly  out  of  phase.  These  data  are 
consistent  with  the  observations  of  Simpson  (1976)  regarding  fingers 
of  high  velocity  fluid  as  mentioned  above. 

3.7  Int.ermi ttency  measurements 

The  intermi ttency  factor  Y,  or  the  fraction  of  time  that  the 
flow  is  turbulent  at  a given  spatial  location,  was  determined  for 
flow  C using  smoke  in  the  boundary  layer  as  the  marker  of  turbulent 
fluid  and  scattered  light  from  this  smoke  at  a given  location  as  the 
detected  signal,  using  the  optical  apparatus  described  in  section  2 
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above.  Due  to  the  large-scale  mixing  in  the  turbulent  region  the 
smoke  is  rapidly  diffused  while  smoke  is  carried  across  the  turbulent/ 
non-turbulent  interface  by  the  much  slower  process  of  molecular  dif- 
fusion. Thus,  the  smoke  was  assumed  to  be  effectively  confined  to 
the  turbulent  regions  of  the  flow.  Fiedler  and  Head  (1966)  verified 
that  values  obtained  by  this  method  agree  with  those  obtained  by 
analysis  of  hot-wire  signals. 

The  signal  passed  from  the  photomultiplier  tube  through  a trig- 
ger circuit  which  produced  positive  rectangular  pulses  when  the  input 
was  above  a preset  discrimination  level.  The  discrimination  level 
was  set  at  approximately  10%  of  the  peak  PM  tube  sign  amplitude. 

The  counter-timer  determined  the  fraction  of  time  the  positive  pulses 
were  present,  which  is  the  intermittency.  Due  to  the  steep  slope  of 
the  PM  output  pul ses , y was  not  very  sensitive  to  slightly  different 
discriminator  levels.  Preliminary  data  taken  in  a zero  pressure 
gradient  flow  were  found  to  be  in  close  agreement  with  the  data  of 
Klebanoff  (1954). 


The  intermi ttency  factor  y of  this  turbulent/non-turbulent 
interface  has  been  previously  found  to  be  well  represented  by  the 
integral  of  the  normal  distribution  curve 

.2  ‘ 


y = — , exp 
/2v  a ■'£, 


1C 
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(16) 


where  5 = y - Y,  Y is  the  mean  distance  from  the  wall  to  the  inter- 
face where  y = 0.5,  and  a is  the  deviation  from  the  mean.  If  this 
interface  is  viewed  as  a wavy  pattern  moving  approximately  at  the 


; j 
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free-stream  velocity,  then  a characterizes  the  amplitude  of  this 
pattern,  y was  in  good  agreement  with  equation  (16)  for  flow  C. 

Table  5 presents  the  characterizing  parameters  which  describe  these 
intermi ttency  distributions. 

The  results  of  Fiedler  and  Head  (1966)  showed  that  the  para- 
meters Y/5*  and  a/6*  are  strongly  dependent  on  the  shape  factor  H 
for  values  of  H less  than  about  2.  This  dependence  is  quite  strong 
as  the  lower  limiting  value  of  H is  approached.  The  minimum  H found 
for  flow  C was  about  1.26  at  about  165  inches.  Upstream  of  the  mini- 
mum in  H for  flow  C,at  130.6  and  149.3  inches,  values  of  these  para- 
meters were  well  within  the  scatter  of  the  Fiedler  and  Head  results 
for  their  "normally  developing"  boundary  layers.  As  pointed  out  in 
section  3.2  above,  no  entrainment  of  non- turbulent  fluid  takes  place 
after  the  minimum  H value  occurs  so  it  is  not  unexpected  that  V/5* 
and  a/6*  downstream  do  not  follow  the  results  of  Fiedler  and  Head. 

The  celerity  and  mean  velocity  profile  measurements  for  flows 
A and  B reveal  that  when  there  is  no  entrainment,  the  entire  inter- 
mittent region  is  moving  with  the  free-stream  velocity.  This  means 
that  there  cannot  be  any  relative  streamwise  motion  of  the  inter- 
mittent region  nor  any  engulfment  of  non- turbulent  fluid  by  turbu- 
lent bulyes.  U (V  - 6*)  and  IJ^a  remain  approximately  constant  down- 
stream. As  for  flows  A and  B,  the  intermittent  region  plays  no  role 
in  the  momentum  transport  after  entrainment. 

Measurements  of  the  frequency  with  which  turbulent  bulges  pass 
a fixed  point  were  made  by  counting  the  number  of  pulses  per  unit 
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Station 

(inches) 

inches 

a, 

inches 

.ZJ 

-< 

X 

o 

1 

aU  , 

xio 

V 

"i  ’ 
1 Y 

Hz 

130.6 

3.00 

0.491 

1.506 

2.465 

47.7 

149.3 

2.55 

0.417 

1.540 

2.518 

45.4 

164.2 

1.93 

0.440 

1.605 

3.660 

50.2 

170.8 

1.59 

0.437 

1 .620 

4.453 

53.0 

178.3 

1.19 

0.306 

1 .759 

4.523 

54.8 

185.8 

0.704 

0.193 

1.562 

4.281 

58.7 

Table  5.  Intermi ttency  parameters  for  flow  C. 


Station, 

inches 

Y 

165.3 

0.02 

171.8 

0.15 

175.2 

0.55 

177.4 

0.78 

182.0 

0.96 

184.2 

1.0 

186.2 

1.0 

- - — 

Table  6.  Surface  intermi ttency  parameters  for  flow  D. 


time  from  the  schmitt  trigger  output.  A serious  problem  with  this 
method  arises  due  to  the  fact  that  several  short  duration  pulses  may 
occur  as  the  probe  volume  enters  and  leaves  the  turbulent  bulge  thus 
giving  rise  to  frequencies  which  are  perhaps  an  order  of  magnitude 
higher  than  the  actual. 

In  order  to  reduce  the  effects  of  this  problem,  the  signal  from 
the  schmitt  trigger  was  processed  through  a low  pass  filter  which 
attenuated  the  short  duration  pulses  responsible  for  the  higher  fre- 
quency. It  was  found  that  the  resulting  frequency  obtained  was  quite 
dependent  on  the  filter  setting,  so  the  filter  setting  was  obtained 
by  visual  comparison  of  the  unfiltered  PM  tube  signal  with  the  fil- 
tered schmitt  trigger  signal  on  a dual  trace  storage  oscilloscope. 
Typical  signal  and  oscilloscope  traces  from  this  method  are  given  by 
Strickland  and  Simpson  (1973).  The  filter  was  set  such  that  there 
was  a single  zero  crossing  from  negative  to  positive  for  each  "sig- 
nificant" peak  on  the  unfiltered  PM  tube  signal.  The  filter  setting 
chosen  for  all  the  results  here  was  300  Hz. 

Figure  33  of  Simpson  and  Wallace  (1975)  shows  a typical  bell- 
shaped  frequency  distribution  such  as  obtained  for  flow  C.  The  fre- 
quency distributions  across  the  intermittent  region  for  flow  C are 
rather  similar  when  the  frequency  is  normalized  on  the  frequency  at 
Y,  n.  . As  for  flows  A and  B,  the  peak  frequency  for  each  profile 

’r 

occurred  at  about  the  location  where  y = 0.6.  The  n.  for  flow  C 

nY 

was  about  twice  the  bursting  frequency  n^  when  tables  3 and  5 are 
compared.  Since  the  technique  used  to  obtain  the  intermittency 
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frequency  requires  a certain  amount  of  arbitrary  judgment,  as  do 
techniques  using  hot-wires,  it  cannot  be  determined  if  the  absolute 
values  of  these  measurements  are  accurate.  However,  since  for  each 
flow  n.  was  proportional  to  n.  , it  appears  that  the  large  motion 

D 1 y 

influencing  the  intermittency  also  influences  the  wall  bursting  fre- 
quency even  downstream  of  the  cessation  of  entrainment.  Strickland 
and  Simpson  (1973)  also  found  n.  to  be  proportional  to  n.  for  an 

D 1 y 

adverse  pressure  gradient  turbulent  boundary  layer.  Thus,  apparently 
this  proportionality  is  present  in  turbulent  boundary  layers  over  a 
wide  range  of  pressure  gradient  conditions. 

For  flow  0 the  intermittency  at  the  surface  was  determined  for 
the  retransition  region  downstream  from  the  flush-surface  hot-film 
signals.  The  taped  signals  were  examined  for  high  frequency  content 
since  considerable  mixing  and  high  frequency  oscillations  are  asso- 
ciated with  turbulent  fluid.  These  signals  were  passed  through  the 
model  330B  filter  with  a 200  Hz-2KHz  bandpass  and  then  through  the 
multiplier  used  as  a squarer  in  order  to  rectify  the  negative  por- 
tions. This  signal  was  low-pass  filtered  at  300  Hz  using  the  model 
3202  filter  and  passed  through  a schmitt  trigger.  The  resultant 
pulses  were  visually  compared  with  the  original  signals  on  a storage 
oscilloscope  for  validity  as  regions  of  high  frequency  content. 

These  pulse  periods  were  accumulatively  counted  for  50  sec.  by  the 
digital  timer-counter . The  fraction  of  time  pulses  were  present  was 
the  intermittency  y.  These  results  are  presented  in  table  6.  These 
same  signals  were  processed  by  using  the  schmitt  trigger  on  the 
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output  signal  of  the  model  330B  filter.  The  resulting  fraction  of 
time  positive  pulses  were  present  was  doubled  to  get  the  intermit- 
tency  since  contributions  from  negative  pulses  had  not  been  recti- 
fied. These  latter  results  were  within  5%  of  the  results  of  table  6. 

The  results  in  table  6 follow  a curve  of  the  form  of  the  inte- 
gral of  the  normal  distribution  curve,  y increases  slowly  from  about 
165.3  inches  to  171.8  inches,  increases  more  rapidly  up  to  177.4 
inches,  and  thereafter  increases  slowly  to  unity  at  about  184.2 
inches.  Upstream  of  165.3  inches  the  flow  had  low  frequency  oscil- 
lations as  shown  in  figures  28(a)  and  29.  It  is  clear  in  figures 
30  that  higher  amplitude  pulses  were  occasionally  present  at  165.3 
inches.  Figures  28(b)  and  (c)  show  the  relatively  high  frequency 
oscillations  in  the  high  velocity  pulses  at  171.5  and  175.0  inches. 
The  clearest  oscilloscope  traces  are  shown  in  figures  31  for  182.0 
inches.  The  higher  velocity  regions  contain  high  frequency  oscilla- 
tions. Regions  of  laminar-like  flow  follow  each  high  velocity  pulse. 
These  velocity  signals  slowly  decay,  not  because  of  limited  sensor 
response  but  because  of  temporarily  high  velocity  laminar  flow  that 
follows  the  turbulent  pulses. 

It  is  clear  that  in  the  low  velocity  wall  regions  downstream 
of  the  cessation  of  entrainment  in  flow  D,  the  flow  became  progres- 
sively more  laminar-1 i ke.  The  magnitude  of  the  oscillations  in  these 
regions  progressively  decreased.  However,  beginning  at  165.3  inches, 
high  velocity  fluid  moved  in  toward  the  wall.  The  average  frequency 
of  occurrence  of  high  velocity  pulses  in  this  downstream  region  was 
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about  25  Hz,  as  determined  by  counting  these  pulses.  Further  anal- 
ysis of  these  data  is  planned. 

4.  DISCUSSION 

Many  aspects  of  flows  C and  D are  similar  to  previous  results 
reported  by  other  investigators  and  by  Simpson  and  Wallace  (1975). 

In  flow  D and  in  the  flow  of  Blackwelder  and  Kovasznay  (1972),  K 
increased  to  a value  which  exceeded  3.8  x 10  ^ and  decreased  down- 
stream to  zerr  . In  both  flows,  a local  maximum  in  C^/2  occurred  at 
the  same  streamwise  position  as  the  K maximum;  the  minimum  in  C^/2 
occurred  downstream  of  the  minimum  of  the  intermi ttency  at  the  wall; 
the  maximum  H occurred  upstream  of  the  C^/2  minimum  and  approximately 
at  the  minimum  intermi ttency  at  the  wall. 

In  all  four  flows  measured  in  this  research  program,  the  behav- 
ior of  the  turbulent  boundary  layer  after  the  beginning  of  strong 
acceleration  was  very  similar.  An  initially  normal  low-pressure- 
gradient  boundary  layer  first  became  distorted.  The  wake-like  tail 
of  the  mean  velocity  profile  became  progressively  weaker  until  at 
the  cessation  of  entrainment  none  remained.  It  appears  that  the 
shape  factor  H is  about  a minimum  at  this  location.  The  bursting 
frequency  and  intermittent  bulge  passage  frequency  downstream  either 
increased  or  remained  about  constant  along  this  length.  The  surface 
spanwise  structure  was  not  appreciably  different  from  a low  Reynolds 
number  low-pressure-gradient  boundary  layer.  Downstream  the  entire 
intermittent  turbulent/non-turbul ent  interface  was  outside  the 
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boundary  layer  thickness  6 and  the  celerity  of  this  region  was  the 
free-stream  velocity. 

Simpson  and  Wallace  (1975)  showed  that  normal  stresses  produc- 
tion was  important  in  the  mean  turbulence  kinetic  energy  balance 
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Here  the  terms  are  from  left  to  right:  advection,  turbulent  diffu- 

sion, dissipation,  shear  production,  and  normal  stresses  production. 
This  last  term  is  of  opposite  sign  to  the  shear  production  term  in 
strongly  accelerated  flows.  The  net  turbulence  energy  production 
is  less  than  the  shear  production  by  the  factor 


F = i _ (u2  - ?)9U/9x 
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(-uv)3U/3y 

which  represents  total  production  to  shear  production. 

They  also  found  that  the  normal  stresses  term  reduces  the  dif- 
fusion of  turbulence  kinetic  energy  into  the  freestream,  or  the 
entrainment  of  non- turbulent  fluid  into  the  boundary  layer.  The 
velocity  difference  between  turbulent  fluid  in  the  intermittent 
region  and  the  irrotational  freestream  fluid  approaches  zero  with 
strong  acceleration.  Thus,  velocity  profile  instabilities  which 
produce  the  eruption  and  eventual  engulfment  of  free-stream  fluid 
are  eliminated.  There  could  be  no  creation  or  merging  of  adjacent 
large  eddies  that  accompanies  the  entrainment  process.  This  also 
explains  why  the  intermi ttency  and  bursting  frequencies  for  all  four 
flows  remained  about  constant  after  cessation  of  entrainment. 
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Simpson  and  Wallace  found  that 

\ - 09) 

described  the  entrainment  velocity  upstream  of  cessation  for  flows 
A and  B.  This  scale  characterizes  the  large  eddy  structure  that  is 
not  only  largely  responsible  for  the  various  turbulence  intensity 
levels  in  the  turbulent  fluid,  but  is  directly  responsible  for  the 
entrainment.  When  the  negative  normal  stresses  turbulence  energy 
production  term  arises  under  strong  acceleration,  the  available  tur- 
bulence kinetic  energy  in  the  outer  region  is  reduced  so  less  large 
eddy  energy  is  available  for  entrainment.  This  equation  reduces  to 
Bradshaw's  relation  for  low  acceleration  boundary  layers  when  = 1 . 

is  given  by  equation  (18)  with  all  quantities  evaluated  at  the 

— 2 

boundary  layer  thickness  6.  Using  the  relations  -uv  t 0.15  q , 

~~2  1 ~2  ~~2  ~2 

u * 2 q » and  v ~ 0.2q  > which  apply  to  entraining  flows,  they 


obtained 
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2K 
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Figure  36  shows  this  quantity  along  flows  C and  D.  Note  that  nega- 
tive values  of  F^  occur  after  engulfment  of  non-turbul ent  fluid 
ceases.  This  does  not  mean  that  turbulent  fluid  is  instantaneously 
reverting  to  a laminar  state,  but  simply  that  progressively  more 
turbulent  fluid  lies  outside  the  6 location. 

Figures  7 and  8 show  the  results  obtained  by  integrating  equa- 


tion (19) 


,(<5  - 6*) 
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For  these  computations  the  relation  (_L|v)rax  ~ 0-3(u  )[;lax  was  used 
"2 

with  experimental  (u  )max  values.  The  results  are  iri  good  agreement 
with  values  determined  directly  from  the  mean  velocity  profiles, 
especially  for  flow  C.  Thus  the  growth  and  decline  of  the  mean  flow 
boundary  layer  is  fairly  well  predicted  by  this  model,  at  least  to 

just  downstream  of  the  cessation  of  entrainment. 

— 2 

The  observed  flat  nF(n)  region  for  u spectra  emerges  after 
some  length  of  acceleration.  The  streamwise  location  where  this 
behavior  is  first  noticed  does  not  seem  to  be  directly  related  to 
the  location  where  entrainment  ceases.  For  flow  A with  K = 2.17x10 
the  nF(n)  spectra  examined  by  Simpson  and  Wallace  had  no  flat  region, 
even  downstream  of  the  cessation  of  entrainment.  For  flow  B with 
K = 3.19  x 10  k,  no  flat  region  was  present  in  the  spectra  obtained 
upstream  of  the  minimum  H,  while  such  a region  was  observed  for 
spectra  downstream.  In  flow  C presented  here,  the  flat  region  in 
some  spectra  began  at  about  137  inches  while  the  cessation  of  entrain- 
ment did  not  occur  until  about  165  inches.  Similarly  for  flow  D 
entrainment  ceased  at  about  132  inches  while  the  flat  nF(n)  region 
began  at  about  123  inches.  A preliminary  exploration  of  this  behav- 
ior is  that  with  strong  acceleration,  progressively  weaker  new  low 
frequency  oscillations  are  produced  by  the  large-scaled  motion  while 
the  higher  frequency  oscillations  are  produced  by  the  breakdown  of 
the  more  intense  upstream  large-scaled  motions.  Further  analysis  of 
these  data  is  needed. 
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In  flows  A and  B,  asymptotic  similarity  flows  were  approached 
near  the  downstream  end  of  the  test  section.  This  means  that  77 /u 

oo 

vs.  n similarity  was  present  for  K < 3.19  x 10  . For  flow  C when 
K > 3.6  x 10  /U^  profiles  decayed.  This  means  that  had  this 

K level  been  maintained  downstream,  true  relaminarization  would  have 
eventually  resulted,  as  suggested  by  Kline  et_al_.  (1967).  Since 
there  is  a much  shorter  distance  from  the  beginning  of  acceleration 
to  the  maximum  K in  flow  D,  the  intensity  n/?  /U^  did  not  begin  to 
decay  near  the  wall  until  about  157  inches.  At  161.6  inches,  as 
shown  in  figure  20,  a considerable  amount  of  the  intensity  was  due 
to  unsteadiness  rather  than  turbulence.  Since  K decreased  downstream 
of  144  inches,  retransition  to  a low  acceleration  turbulent  boundary 
layer  began  at  about  165  inches.  In  this  region  the  flow  possessed 
a highly  laminar-like  behavior,  so  -uv  was  extremely  small.  Evidently 
77  decays  considerably.  As  mentioned  above,  the  surface  spanwise 
structure  upstream  of  the  onset  of  decay  of  77  /U^  was  not  appre- 
ciably different  from  that  of  a low  Reynolds  number  low  pressure  grad- 


ient boundary  layer.  After  the  onset  of  decay  of 


77 


/U  , X’ 


increased.  Evidently  the  upstream  values  in  these  non-equilibrium 
flows  persist  downstream,  so  when  normalized  with  a higher  U^/v,  X* 
is  higher. 

In  the  relaminarizing  nozzle-type  flow  of  Schraub  and  Kline 
(Kline  et  al . 1967),  A*  was  reported  for  the  region  of  maximum  K to 
be  of  the  order  of  200  for  K = 2.75  x 10  ^ and  of  the  order  of  300 
for  K - 3.25  x 10'6.  On  the  surface,  it  would  appear  that  they  found 
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A+  to  vary  directly  with  K.  However,  their  estimates  of  C^/2  for 
these  two  streamwise  locations  seem  high  by  about  a factor  of  two 
when  compared  to  the  C.^/2  behavior  in  flow  D and  in  the  Blackwelder 
and  Kovasznay  (1972)  nozzle  flow.  In  that  case,  values  of  about  150 
and  220  for  A*  would  result,  respectively,  and  would  agree  reason- 
ably well  with  the  flow  D results  in  the  vicinity  of  the  maximum  K. 

Schraub  and  Kline  did  not  present  Az  results  downstream  of  the 
K maximum,  so  the  data  from  flow  D are  the  only  available  on  this 
downstream  structure.  A*  is  still  of  the  order  of  300  while  K 
decreases  below  2 x 10  After  retransition  to  a low-accelerated 
low-pressure-gradient  turbulent  boundary  layer  begins,  A*  is  an  order 
of  magnitude  greater.  This  behavior  is  as  yet  incompletely  explained, 
so  further  analysis  of  these  data  is  needed. 

5 . CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  can  be  made  about  turbulent  boundary 
layers  subjected  to  strong  acceleration. 

a.  The  parameter  K must  be  greater  than  about  3.6  x 10  6 for 
the  decay  of  the  upstream  generated  turbulence. 

b.  For  practical  use  in  a nozzle  only  one  short  region  down- 
stream of  the  maximum  K will  possess  a laminar-like  behavior  with 
significantly  lower  Cf/2  and  heat  transfer  coefficients. 

c.  The  large  eddy  structure  of  the  outer  region  governs  the 
bursting  frequency,  the  intermittent  bulge  passage  frequency,  and 
influences  the  wall  flow  behavior  downstream.  These  frequencies 
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seen  to  approach  constant  values  after  the  cessation  of  entrainment 
of  non-turbulent  fluid. 

d.  The  wall  spanwise  structure  appears  to  lag  behind  local 
conditions  and  to  reflect  upstream  flow  behavior.  A*  seems  to 
steadily  increase  to  about  300  even  after  K has  dropped  below  its 
maximum  value  in  a nozzle  flow.  After  retransition  to  a low  accel- 
erated turbulent  boundary  layer,  a large  A*  of  the  order  of  3000  is 
observed. 

e.  The  proposed  modified  entrainment  model  accounts  for  the 
reduction  of  available  turbulence  kinetic  energy  by  the  negative 
normal  stresses  production  term  of  the  turbulence  kinetic  energy 
equation. 

~2 

f.  A nF(n)  spectral  distribution  of  u for  n > 1000  appears 
to  possess  a flat  region  for  laminar  like  boundary  layers  at  large 
K values. 

As  a result  of  the  large  quantity  of  experimental  data  provided 
by  the  research  program,  the  following  recommendations  for  future 
work  are  suggested. 

a.  Further  analysis  of  these  data,  some  of  which  was  not  pre- 
sented in  this  report  because  of  the  lack  of  time.  Afterwards  a 
manuscript  should  be  submitted  for  journal  publication.  At  that 
time  a more  intelligent  assessment  for  future  work  can  be  made. 

b.  A generalized  computational  effort  incorporating  flow  models 
reflected  by  all  available  data  should  be  undertaken.  Communication 
between  the  senior  author  and  predictors  has  begun,  but  no  results 
are  yet  available. 
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APPENDIX  C 


ERRATA  FOR  REPORT  SMU-l-PU 


1 ine 

bottom  "96  inches"  should  "93  inches" 

6 "3.1672  x 10'4(x)"  should  be  "3.1627  x 10“4(x) 

6 "6.5931  x 10‘6(x2)"  should  be  "6.5931  x 10  7(x2) 

"1.1944  x 10"5(x2)"  should  be  "1.1994  x 10  5(x  ) 


8 


AD-A 


A037  Ml 


UNCLASSIFIED 


PURDUE  UNI V LAFAYETTE  IND  PROJECT  SQUID  HEADQUARTERS  F/6  20/4 

LAMINAR I SCENT  TURBULENT  BOUNDARY  LAYERS:  EXPERIMENTS  ON  NOZZLE  —ETC(U) 
JAN  77  R L SIMPSON*  C R SHACKLETON  N00014-75-C-1143 


1 


SQUID-SMU-2-PU 


NL 


END 

DATE 

FILMED 

4-77 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURLAU  Of  STANDARD^  I'tfH  A 


I I 


E i 


■ *0 


>>C\J 
4-  CM 

o 


03  03  - 


4-  03  C -r- 


03  4~»  4-  — J 


_ E 
00  4- 

r-  o 

O 4- 

c 


i > LD  c 
3 •—  O O 
-QCVJU 


O C — 

X •«“  03 
^ Oiu 
_J  C- 
■r-  C 

in  > X 


U 03  ' 


^ c u c 


03  3 - X 


03  X 

■O  o 
o ce: 


o - 03 
a ct- 
<o  oj 
<U  4- 
3 cr  •• 


i.  ao  ai  o 


: ro  q: 


a»  c cnh-  c_> 


. co 


<—  S-  2 3- 
<o  o ai  h- 
>u.z< 


a> 

3 co  ■*-> 
c in  4- 
=)  o •• 
qj  xz 
— T3  3 3- 
oj  o oj  t— 
> o z:  < 


03  4-  i- 
c 4—  a> 
QJ  O 3 
in 
4- 
•a: 


i_>  +->  4- 

o»  h—  • o 
in  •«-  o 
o • >> 

m Q 4- 

■*->  -o 


? ° 

a;  <u 


C C - 4- 


>^•0 

c c - 

C 03  4- 


4->  E > ^— 
03  O O 3— 

u o o «t 


03  03  C X 

E O •*- 

3 CD  4->  _J 
U C73  C7> 

O C C 
CO  03  -i-  • • 
X X Z 

0) 


p X «3h- 


03  UJ  3 C 


03 


X CO  >>cm 

<->  r--  > 

<3-  03 


X 

> CM  O 
OJ  CM  O 
2=  CM  <_> 


>>0  E 
i - o>  i- 
o CO  o 

4->  O 4- 
03  CM  C 
4-  — ■ 


CT> 


■a 


au  p oi 


r—  +->  03  C 


03  OO  tn  03 

on  E 
4-  oj  E 
03  2T  O 


03  X O > 


O 4-  Q O 


-O  • <— 
03  O 03 
— I • U 

Q -r- 


Q 4- 


O 03  *J  C 


03  O C 

03  r—  ♦->  2 
CO  f-  C3 
03  f~-  2 


03  C 


u in  i/i  i — 


4-  «3-  co  c 

4- 

o 


f-Q.Q< 


I X in  H- 
■ <V  03  3— 
. C3  3 <C 


03  03  3— 
>3< 


O C 

->x  x 
» co  u 

J O 03 


>>X  03 
> CO  i_ 
03  O XI 
znj-r 


vHJur 
E • o O 
03  4-  O 4-  *— 
-4-3  O 03  CO 


• o 

> C_3  •»- 

• i/l  O 

■ Q r-  ro 
) 3 co 


C CO 


>,-*->  C on  X 

x c o oj  «— 

03  -M  X <t 


cn 


■ *->  c 


u i 


co  cn 

o 4-  — 
<->  CL  et 


. 4-  *-  •• 
oir  z 
• X in  3 — 


*—  X 
03  03 
> Q 


X JZ 
on  h- 
03  3— 
3 < 


>,+-»  c u or 

X C O 03  — 
03  ■*-*  t — e£ 
c E 03 
c 

< i-  . . 

<o  X z 
■ — x in  t — 

03  03  03  ►— 

• Q 3 C 


3 O 03  CO 

: cn  •<- 

3 r-  -O  . 


13  X CJ>»- 


C i-  JD 
03  -*  a.  -r- 
C 03  l/l  -I 

-a  o 


03 


3 =5  3-  3- 
: co  u < 


o «a- 


>>0  CM 
s-  a>  cd 

O CO  CM 
4->  O 
03  CM  03 


CM  c 


T3  T3 

: c o 

3 03  O 


03 


u »•  c 

C C X 
03  O •«- 

03  ■*->  I 

on  cn 


l/)  O.  Q.  •• 


03  03  ; 


o o 
CO  >,x> 


in  . 


03 


£ *-> 
03  S~ 

^ c 


cr  x o Q 
o 

E CQ  03 


ffi 


S- 


OO 


C C H- 
03  C C f— 
><<< 

5 


o 


oo  • co 

4-  Q co 
03  O O 
O 

C *J  C Q 
03  C O CC 
C 03*30 
T3  E Oi 
M 4->  C 
O ^ -r-  •• 

oj  r z 

<—  o.  «/» y- 

03  03  03  I— 
> Q 3 < 


03 


a>  c * v. 

4-3  03  >,C3 
cn  E 03 

o +-»  s- 


•—  a.  c h- 


03  03  Ol- 


03 


a j 


lD  O 
S-  LD  •>- 
03  LO  cn 
•a  co 
c cn  3 


03 


OJ  >>vx  03 

§>  CO  i. 
03  O X 


cn  x • 
E*3u  oi 
03  • a 

4->  4-  Q -r- 

cn  o c 

>>  'X 
40  4->  C U 
c o 03 
Q.  03  4-3 

— E cn 


03  03  03  cn 

uu  ua- 


<—  x tn  h- 

03  CJ  T3  H 
> O ^ C 


o 

- <*■ 

</)  cn 

ooi  i/i 

— - cn  x 
0*->  3 

0 3 Li_ 
X 03  ro 
U C -r-  . 

40  O C LU 
i-  4- 

03  03  O C 
4-3  <C  4-  03 


cn 


<—  c 
CO  o 

LTl  -f- 


<_>  -r- 


o 

>»  to 

03  -a 

cn  •— 
"O  s-  03 
03  03  C 

o o o 


a cm  dj  u 

CMU-r 
-O  CM  C 
--  CO  - X 

3 r— 


X O X 

o -*-> 

E CQ  03 
03  U •• 

x • c z: 

O O I — 

o • s-  *— 

S-  Q_  Q_  <x; 


03 


X C t— 


03  C C C i- 


•'-OO  U »— 


03 


C cn  *—  /—  m 
03  •<-  03  r—  o 

O > O -r-  x 


cn  Q 


03 


_ >,-x  S-  o 

Oi-  (CQU 
03  4->  _J 

03  cn 

3:  •<-  os  • • 

E c z 

— 03  -r-  t— 


03  - 


>uu< 


> u u < 


X X 3 

03  O U 

4-  *—  4- 

fO  *—03 


<n  : 


4-  03 

- cu  •<- 

in  4-3  c 

4->  in  os 

U (U  > 

03  X r- 

■r)U  >, 


in 


4-  -M 

x in  c 

O 03  03 

4-  3 O- 
03 


O 4-3 

X 03  C 
u — ' "f—  03 
40  *3-  c -O 
CM  4-  C 
03  *—  O 03 

4J  CM  4-  4-J 


3 0)*-^ 
T3  n 1T5  4. 
03  O O 03 
4-0  X 
cn-—  - 3 

4-3  >>40 

cn  >>  03 
0 4-4- 
Q.  03  4->  • • 

4-  03  Z 

— x c h- 

o---  OH 


3 *— 
"O  03 
03  O 


X 4-3 

<J  03  • 

4-  Oj  03  '"O 
03  4-  -r- 
03  4-)  C X 
cn  i/>  4-  X 
0)  O <— 
C£L  C 4-  O 
03  ■<-  T3 
•“  03  *—  3 
US  4-  03  CC 
> O O 


4-3  X 
03  O 
03  X 


03  4- 


03  ■ 


03 


i/i  in  i/i 

03  -r-  4- 

a j;  o u 
4-  CO 
*—  03  •*-  C 
03  ■ — ■ — 03 

>C_3r-  P 


03  — 


*—  C f— 
03  O 3— 

> 51  < 


cn  03  Q 
4-  03  C 
C?  LU  03 

■o  •• 

CQ  UZ 

o'  ro  m »— 

■<-0  <0  3- 

4-  r-  (1< 


Z X 


U X •*“  »— 
-r-  ro  X 3— 
4-  X O 


I • 


'IMP»  ... 


wmmmm 

* 


* 


I 


O CT*  Qj 

c na 

X r- 

O C\J  r— 

vU  O 03 


>,  atco  - a. 
4-»  E ^ O 
■r  *j  M cn 
i/i  n> 

c.  to  os: 

0)  Q.  X ••—  I— 

> qj  o x t— 

•r  Q CO  U < 


cno 
c co  s- 

CTiLU  <X>  qj 

c O -X 

tO  t-  to  u 
•r-  s-  a*  ra 

o a»  i-  m 

C D OJ  </> 

C C •«-  • 

*—  -r-  UJ  O CO 
— CT  C 
C 4~  • 

UJ  O I — o 
<4-  ■— 

O 4-  +j  ~-*  • 

o c s- 

>,  U - O 
■*->  QJ  E O 

•»-  C7*  4->  CD 

1/1  0)  i-  *0  •• 

S.*-  uz 

QJ  <—  Q.-r-t— 
> O O X 
•*-  C_>  CL  t_J  <C 


co  13  cn 

i o c.<c 

I C\j  o 

l-o  s-  c 

I fo  D_  O 

0-0  -r- 

) IX J *tj  aj 

• to  •—  u E 
> C S- 


3 U UW 
■*->  QJ  3 S- 

c s:  ♦-  a> 

OJ  C X 

x:  o-  a,  c 

c v o: 


•*“  ■*->  C7> 
to  J-  C •• 
S-  *T3  r-  Z 
QJ  Q.  X t— 

> qj  a,  t— 


'f-u  m i. 
O +->  <TJ 
*—  *->  2 
0*0  0)0 
+->  o to  z: 


to  O to  J* 
c to  o 
•—  4-  «o 

o s:  -o 


■*->  c a>  i_ 
C d)  CiO 
to  E T3 
3 4->  *r- 

x s-  i.  •• 

U (DUZ 

*o  a Et- 
to  qj  aj  h- 
to  O O c 


<\j  -m  x_ 

Of  I CD  O 

*->  LU  lO  Q. 
3 3 QJ 

+->  «3 - X QC 

■r-  r—  u 

+->  fO  «— 

(A  E to  1X3 

c o to  o 

•—  O *TJ  •*- 

a:  2:  c 


+J  (/)  (U  Q; 
O)  CD  Oil — 
to  •«—  "O 
3 i-  -r- 
r n3  S-  •• 
u 1.  nz 
T3  -O  E ) — 


00  CO  •*- 

4->  o * O! 

'J  (U  c 
O)  I OILU 
to  o *o 


U E -O  ^ 
*T3  C E t— 
to  o *T3  K- 
«i  a u < 


>,  OJ  o 

4JU  0£ 
•<-  -0-0 
oo  x c 
i.  U J-  UJ 
QJ  i-  O C/0 

> *T3  r— 

■f  qj  o • 

C U)U< 


*0  c -r-  . 

+->  -r-  r—  fc. 

co  S---2: 
qj  o 
o qj  C_> 

-o  c 

<TJ  •*-  4->  Z 

i-  oi  i-  J— 

o c o \— 

r—  UJ  U.  C 


>»•*- 
03  X 
O 0) 

*—  <D 

o c 

c 

x oicvj 
u c ro  C 

(Dune 


ZD  ' — i-  O 
O 0)  >- 
s_  o c.  ; 


>,x  -x  s- 
4->  o i-  o 
•f—  O to 
to  4-  >-  to 
i-  O QJ 
QJ  2 4- 

> 4->  O)  o 

•<-  C z S_ 
C QJ  O- 
D E ' 

+->  <v 

I—  i-  u •• 
>-  UJ  ITJ  Z 
QJ  Q-x:  h" 
c aj  -m  |_ 
i-  a •—  c 
o 


>,**- 

cn  s-  c 

O QJ  *0 

•—  0)  i- 

o c o 

C -r-  cr>  U 
x:  cno  x 

u c o 


X o 

U >1  *— 

QJ  s-  *— 

►—  o CT* 


a>  E C 3 

4->  QJ  C- 
3X0  1- 
-*->0  4-0 

•r-  -r-  to 

4-J  4-  r—  to 
to  O ^ QJ 
C 0 4- 

— • 4->  o 

C " 4- 


S-  1-  -O  •• 
o <xj  n z 

4-  Q.  «/l  ►— 

•r-  qj  *o  h- 

r-  Q Q < 


o i-  a)  flj 

o > •*- 

Ql  X ■*-  C 

■P  ^ c c 

3 _J  o O 

4- >  4- 

•r-  C QJ  •*-  >, 

■*-»o>*—s- 
to  O *TJ  fO 

c to  S-  O i- 

— O X 

3 »•  -t- 

*T3  Cl.  -X  *T3  _J 

o *T>  c 

C C.  O QJ 

5-  CX.  T3  • • 

O O *tJ  z 

4-  4->  o to  H- 

•f-  QJ  CO  <T3t— 

r-0*t  a< 


cn 

vo 


O 


c\j 

r>* 


ro 

r-*  f'** 


un 


vo 

r^. 


CO 

r^. 


•<-  3 — 


•r-  03  03 


!£. 


s i 


to  o 

vo  to 
03  to 
51  QJ 
«•- 
» *•  O 

QJ  t- 
I 0>£X 


!«^ 

> 03  ►— 

: c_>  < 


U •*J  E 
-0  0.0 
00  O)  o 

oO  Q QC 


3 <u 

*->  -f-  c 

•r-  -O  S_  r- 
■M  t- 

tO  03  w 

C > tO 

— • QJ  to  • 
*—  •-  < 
x:  3 z 
o o • 

s-  CO  - H— 
03  >> 

QJ  t-  4-»  . 

0O  <U  r-  1_ 

(U^UQ 

QC  r- 

O 00 
» 2» 


■f-  > *—  x : 


2 1 


4-i  T3  • O 
■r-  03  >,U_ 
4->  O L. 

wa  3 • 

c -O  s- 

o 

Qj  00 
-*  *—  QJ 


s-e 


Qji-.cc 
£(/)  HJ4J  OJ 


o •—  o — • 

03  -r-  i_ 

>,000X1 


QJ  03  OO 

S-  o 

4->  i.  "O  *-> 

O TD  ns  C C 

Z V V. 

o X3  u 
*♦-  O — ' Z 

o a)  _j  a> 


00  qj 

03  QC  to 

E 4-> 
vo  o i-  u 
in  u o cd 
on  u **- 


qj  oo  vo  ' 

E c ^r 

03  •*- 

Q V (D  | — 

qj  c 

QJ  QJ  03  ■*-> 
l-  C •«“  t- 
4->  •»-  TD  03 
O cn  C 3 


O £ 
>,  03  £ 

+->  a>  o 
•<-  cn 
to  a>  qj 


03  (X 

QJ 

Q r- 


>,LlJ  3 
«-»  C 

•*“  — QJ 
00  03  > 
i-  U C 
<U  •£ 

> E -C 

•r-  <U  4-> 

c .c  a> 

D U r- 

Qj  *4-  4-J 

•*->  o to 

03  QJ 


x cn 
<u  c 
> •<- 
-r-  -O 

c <— 

13  -f- 

3 

qj  ao 


03  QJ 
•*-  TD 
C •«- 
03  VO 
> QJ 
f—  L 
>>CL 


m 4J 

^ x 

QJ 
O -Q 
■r-  O 
-C  QC 
O 


£ .. 


•r-  T3 
C *—  ( 
03  O 


•>  *J  4-  O 
* *r-  O 


■Q  z 

to  vo  Z 

a»  -m  z 

o 

03  Z 

X - 5 Q.  Z 

fc-  T 9 

1-  Z 

1-  i-  z 

to 

+■> 

3 Z 

QJ  E 

E t- 

QJ  LO  C 1— 

X Q.  X to  >— 

>- 

QJ  T3  H- 

QJ  L.  03  03  K— 

QJ  v.  E 

4-»  h- 

QJ  — -»J  I— 

ao 

h~ 

Q.  O 

03  f— 

1CM  OJI- 

aioo3»— 

-C  ■—  >— 

>Q)vrh 

> QJ  QJ 

O h- 

> O O b- 

o 

QJ 

O h- 

O 

o < 

■o  x:  «t 

Jaco-K 

3 O «* 

•r-  Q.  Q U < 

•flrtZ 

z < 

- o Z «* 

Q — 

o c 

qj  ac 

QJ 

QJ 

c 

c 

c 

XT 

XT 

X 

Z 

Z 

33 

13 

3D 

o 

►— 

m 


4->  4- 
CD  03 
U O 


03  <_> 

•r-  «->  c 

cd  oi  - 03 

O Ol  T3  L 

r—  4-  4.  U_ 

O 4-»  o 

Cl/)1*-  • 

-C  -M  4- 

u c i~  t: 

03  •*-  OJ 

►—  03  :n 


aiOcoh- 
*->  O ^ H 
•r  kl«t 


1-  >» 
. <D  ns  X3 
QJ  -O 
C C •*- 


l/l  CO  . 

aioi/) 
. o co 


a>  cd 

•r-  4~>  C 
CD  d)  * 03 
O 03  TD  *— 
— 4-  4.  < 
O -*->  O 
Clrt1*-  • 
.c:  ■*->  4. 

C3  C 4-  O 
a;  -r-  <o 
H flJX 

s: 

TD  wz 
(DO  i/)h 
■*-»  O 03  t— 
■>—  «fU< 


<D  4-  U O • • 
i.  <o  anz 
a;  cxs:  f— 

> OJ  Oh 

•*-  Q —I  C 


O a; 
<—  < 
o 

CD  »4- 


aj  a. 
> a; 
•*-  Q 


- o -a 
u td  c 
oo  03  •<- 


•«-  <u 

CD TD  •• 
C *—  Z 
LU  Z3  ►— 


4- 

03 

a; 

c 


OJ  00  O 

u ^ to 

03  e 

CQ.  03 

in  c TD 

CD  O 03  <C 
•*-  4.  CD 
-C  OJ  • 
U<  ru 

•r-  O 


u ach 
> 0)  Ch 
•*-  Q < <t 


03  U 

CO.  -r- 

03  ID  C CO 
CD  O 03 
•r—  4-  CD  C 
-C 

u < r *j 


0 <— 

CD  o 

>»cm  c 

1 4-  LT3  03 
*->  ic  0J 

•r-  tO  Q 
4.  *f— 

3 E •-  >, 

O 03  4-  C 
lO  . 


•r-  0 4-  •<— 


lO  CD  O . 


03 


>,  CJ  o o 

■M  E -Q 

•r-  4- 

l/l  1-  < '• 


d Q.  c t- 
> d Ch 
Q <t  C 


Z ^ V)  c 


>,  HI  03  Q 
4->  E -r- 

■stci .. 

4-  03  3 Z 
QJ  Q.*—  h~ 
> 03  O h- 
■r-OU  <t 


CD  O •»—  03 

OJ  C 5 

C Q.  4-0 
4-  ID  0 4- 

03  O U-  CO 

-C  4_  •»- 

4->  03  f—  . 

3<Ji  iDV 

o s-  c_> 

to  * 


03 

o a_ 


to  CD 
r-  ^ 1_  C 
id  t d< 

4-  03  > ; 

03  CX  -f-  id  I 

> <u  c o » 

'rODj. 


CD 

dM  4. 
C O CD  QJ 

o •*—  c 

*->  C (D 

O)  03  C i/I 
C -C  O QJ 
•r  U*J(5 
£ d O) 
l/l  Z C . 
03  •*-  CO 

3 ^ -C  • 
O (D  Li. 


>.  aj  - Q 
+j  E a> 

Ifl  t 4J  •• 
4-  03  4->  Z 
Q)  Q.  03  h 
> QJ  QJ  I— 
Q l/)  < 
C 


C >> 

-C  4-1 
U *r- 

<D  to 
-*-»  4-  t 
>»  <U 

•—  > r 

o — 

Q.  C 
ZD  -i 
03 


03  *—  03 
0-0 
dio  4. 
CO  CO 
C\l  - 
CD  O 

c 

•r*  03  • 

4.  -r-  U_ 
< 13  C 
QJ  4- 
C CD  d 
4.  +J 
CD  •!— 

C > 03 


Qj 
C 4-> 
•f  03 
CD  4-> 
4-  on 


Q-r— 

QJ  tO 

Q O 
*3-  to 
CD  CM  Qj 
c tr> 

r 

ai  c 

QJ  -n-  . 


c •*-»  x:  o z 

o 03  h- 

QJ  r-  I— 


c J 

-c:  -*->  •» 
u t 

a>  id  i 
-*->  4.  L 
>,  QJ 
r—  > r- 
o -r-  r 


•r-  QJ  03  -iC 


•r-  QJ 
C to 

o a> 


a> 
o i> 
x>  c 


n o 4-  ►- 

: o o t— 

j U.  U-  < 


ss 


O DnI 

»D  "O  >> 
3 0*— 
Q.  *— 
E *•- 

T)  ID  3 


(D  *—  4- 

— — O 

— O 

■r-  O 


•*-  ID 

4-  <U 
*->  OC 


& 


CO 

CM 


• z 
>» 

CO  T3  • 


^ x 4-i  a 

UJ  03  o O 


4-  4-  O .C  »— 
QJ  O *U3— 
c o a.  ui  «t 


— CO 
o C'J 
c *— 

-C  I 

U to 

QJ  O 

•“  I 

<4-  o ID 

o oc  d ■*-* 

3 4-1 
QJ  * c d 

4-J  4-1  QJ  to  to 

3 C > 3 tD 


4_>  4->  tO  03 

to  4.  4->  lO  C 

c ti  *j  id  r 

— a.  qj  03  o 

0)  ID  z o 

tD  Q 3 
4->  £ • • 

4-1  >,  U OJ  4_ 

Qj  4.  03  CD  Q 
i/)  4->  to  T3 

3 tO  tO  -r- 

u i^az 

03  QJ  E ►— 
tox:r^o3  3— 
^UN  U< 


•r- 

CD  CM 
C CO 
UJ  CO 


tO  C CD  tD 
4-  03  •*“  tO 
dr  r o 


s:  c 


QJ  O CD 
4->  c 
03  «-»  *-  • 

4->  C to  4. 
C Q 


1/3 


03  4- 


: QJ  03  ► 
i a uj  • 


>>  >,  c 

4-1  4-  03 
•<-  O > 


> O C CO 


ZD  I Q.  Q. 


03  4-  4.  CD 
4-»  03  03  "O 
OO  Qj  Cl.  uj 
to 

03  0)  >>  • 

•-  a:  4->  u 

c •—  o 

03  TD  tO 
> QJ  4. 

03  • • 

>>—  > z 

to  CD.  - »— 

c a.  c »— 
C < ZD  *t 
03 


>-  .*-  IT) 


03  UJ  ' — ^ 

03  tO  4. 

t*_  O O 03 

O 03  -r-  ^ -O 

CD.  C 4-  Q3 

d to  03  O — I 

4-»  O ^ >- 

3 4 U *— 

4-1  03  0J  J QJ 

— **  51  OJ  3 

4->  Z E 

tf)  ^ 3 03 

C O 03  ' 1/3 

— •*-  o> 

4-1*—.— 

u c Q.  03  4. 

03  Q.  T3  £Zi 
c E<  CD 

x:  c 

« i u -o  - •• 

E z 

Q.  03  4-  »— 

^ A 


Q> 


QJ  QJ 
> •«- 
u 

C tD) 


■r-  O 


g! 


- LO 


1/3 


U.  < 


03  LO  O 

0.0  4-1 

m to 
. *3-  C 
CD  -C 
O 
O 


O 

C CM  — i 

TD  -C  l/3 
4-1  u r**  4. 
to  3 03  03 

•r-  r—  )—  CD 
T3  U.  to  O 
O *4-  03  OC 
X TD  O X 


iD  ' — ) — 


d c a 

>,  0)  4. 

4-»  C O to 


D OT3Q 
•p-  4. 

TD  C O 
4-  03  ^4-  • • 

o a c z 

•u  a 03  ►— 

c o»  4->  »— 

0Jil/3< 


03  U*> 
0.0 
QJ  CO 
O *3- 


•*“  03  •*“ 

d C3C 
>,d  4. 

*-*  c o . 

«4-  O 

tO  CD  •*- 
4.  C *—  • 

QJ  UJ  *0  1/3 


U-  U O)  )— 


03  LO 
0-0 
QJ  CO 
O *3" 

CD 

CD  U- 
C r— 
•*-  03  03 
4-  - OO 
O)  C 

>,  QJ  4-  >, 
4J  C O 03 
-r-  •*-  U-  C 
to  CD'*—  T3 
4-  C *—  •*- 
QJ  UJ  03  t/3 


O £Z  C Z 
*4-  O 03  H- 
C Qj  4-i  ►— 
03  2E  1/3 


CM 

O 

CD 


03  •— 
•r-  QJ 

E % 

O 03 
^ O 
•*-  4- 
— CO 


U -X  •“ 
lO  4r  . 

03  O O 
CD  O.  03 
C 03  • 

to  •*-  0)  CO  4- 

E 4-  U O 

03  03  03  O 

■*->  03  Q.  TD 

lO  C to  C .. 

»>•«-  O Z 
1/3CDQJTD)— 
C C dh 

3uoa«t 


Unclassi fied 


SECURITY  CLASSIFICATION  OF  this  P ace  (When  Detm  Entered) 


REPORT  DOCUMENTATION  PAGE 


. REPORT  NUMBER 


RF.AD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


SMU-2-PU 


4.  TITLE  (end  Submit) 

LAMINARISCENT  TURBULENT  BOUNDARY  LAYERS: 
EXPERIMENTS  ON  NOZZLE  FLOWS 


5.  TYPE  OF  REPORT  4 PERIOD  COVERED 


Technical 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS  ••  CONTRACT  OR  GRANT  NUMBERS 

Professor  R.  L.  Simpson  and  Mr.  C.  R.  Shackleton  N00014-75-C-1 143 , 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 


Department  of  Civil  and  Mechanical  Engineering  k 

Southern  Methodist  university 

Dallas,  TX  75275  


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Project  SQUID,  Thermal  Science  and  Propulsion 
Center,  Purdue  University 

l.lor  f I 3 f T M /!  7 0fl7 


4.  MONITORING  AGENCY  NAME  4 ADDRESS!'//  dllferent  irom  Controlling  Office) 

Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 


16.  DISTRIBUTION  STATEMENT  (ol  thlm  Report) 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  & WORK  UNIT  NUMBERS 

NR- 098- 038 


12.  REPORT  DATE 

January  1977 


13.  NUMBER  OF  PAGES 

87 


IS.  SECURITY  CLASS,  (of  th/m  report ) 

Unclassi  fied 


IS*.  OECL  ASSIFICATION/OOWNGRAOING 
SCHEOULE 


This  document  has  been  approved  for  public  release  and  sale;  its  distribution  ; 
is  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  the  ebetred  entered  In  Block  20,  II  dlllortnt  Irom  Report) 


l».  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  rereree  tide  II  neceeemry  Mid  Identity  by  block  number) 


Bursting  frequency 
Cross-correlation 
Dissipation 
Entrainment 
Hot-film  anemometer 


Intermi ttency 

Laminari scent  boundary  layers 
Turbulence  intensities 
Velocity  fluctuation  spectra 
Vi  sco 


joAaBS  TRACT  fConitnue  «i  r*mM  mide  If  ntcMaarjr  mod  Identity  by  block  number) 


Because  a turbulent  boundary  layer  in  a nozzle  undergoes  strong  acceleration, 
a laminari scent  boundary  layer  with  the  benefit  of  a lower  surface  heat  transfer 
rate  can  result  for  some  short  distance.  Here^and  in  a previous  report  (Simpson 
and  Wallace,  1975)^evera'l  phenomena  whicTT  accompany  laminarescence  produced  by 
strong  acceleration  are  examined  for  nozzle- flow  and  sink-flow  accelerational  dis- 
tributions, respecti vely. ^Several  different  type  measurements  of  the  structure 
of  two  nozzle- type  flows  are  reported  to  determine  how  an  initially  normal 


DD  1473  EDITION  OF  I NOV  49  IS  OBSOLETE 


Unclassified 

MCUMTV  CLASSIFICATION  OF  THIS  PAGE  fBw  Dmtm  Entert 


SECURITY  CLASSIFICATION  OF  THIS  P AGF.fWian  Data  Enf.r.d) 


20.  turbulent  boundary  layer  approaches  the  laminar-like  state,  including  mean 
velocity  and  Reynolds  stresses  profiles,  spectra,  turbulent/non- turbulent 
interfacial  structure,  and  wall  bursting  and  sublayer  spanwise  spatial  struc- 
ture. 


As  a result  of  these  experiments,  it  appears  that  the  surface  skin- fric- 
tion is  not  reduced  to  laminar  values  in  sink  flows  unless  K(=vUo^2dUoo/dx)  is 
greater  than  about  3.6  x 10'6.  In  nozzle- type  flows,  K must  also  be°°greater 
than  this  value  over  a short  distance  in  order  to  produce  a short,  relamina- 
rized  region  downstream.  The  large-eddy  structure  of  the  outer  region 
governs  the  bursting  frequency,  the  intermittent  bulge  passage  frequency,  and 
influences  the  wall  flow  behavior  downstream.  After  the  cessation  of  entrain 
ment  of  free-stream  fluid,  these  frequencies  approach  constant  values.  The 
wall  spanwise  structure  appears  to  lag  behind  local  conditions  and  to  reflect 
upstream  flow  behavior.  After  retransition  to  a low  acceleration  turbulent 
boundary  layer  downstream,  much  larger  spanwise  scale  structures  are  observed 
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nF(n)  is  constant  for  laminar-like  boundary  layers  at  large  K values. 
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